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ABSTRACT 

We present semi-analytic dynamical models for giant molecular clouds evolving under the influence 
of HII regions launched by newborn star clusters. In contrast to previous work, we neither assume that 
clouds are in virial or energetic equilibrium, nor do we ignore the effects of star formation feedback. 
The clouds, which we treat as spherical, can expand and contract homologously. Photoionization 
drives mass ejection; the recoil of cloud material both stirs turbulent motions and leads to an effective 
confining pressure. The balance between these eflfects and the decay of turbulent motions through 
isothermal shocks determines clouds' dynamical and energetic evolution. We find that for realistic 
values of the rates of turbulent dissipation, photoevaporation, and energy injection by HII regions, the 
massive clouds where most molecular gas in the Galaxy resides live for a few crossing times, in good 
agreement with recent observational estimates that large clouds in local group galaxies survive roughly 
20 — 30 Myr. During this time clouds remain close to equilibrium, with virial parameters of 1 — 3 
and column densities near 10^^ H atoms cm~^, also in agreement with observed cloud properties. 
Over their lives they convert 5 — 10% of their mass into stars, after which point most clouds are 
destroyed when a large HII region unbinds them. In contrast, small clouds like those found in the 
solar neighborhood only survive ~ 1 crossing time before being destroyed. 
Subject headings: H II regions — ISM: clouds — stars: formation 



1. INTRODUCTION 

Giant molecular clouds (GMCs) are the primary reser- 
voirs of molecular gas within spiral galaxies. Star for- 
mation is tightly correlated w ith the molecular column 
density within spiral galaxies ijWong fc B litz 2002), and 
is therefore controlled by the formation and evolution of 
these giant clouds. In this and a subsequent paper we 
develop the theory of GMC dynamics and present semi- 
analytical models for GMC evolution. We rely on simpli- 
fying assumptions about the structure of the cloud and 
the properties of the surrounding interstellar medium in 
order to focus on clouds' global energy budget and dy- 
namical state. 

Observations show that stars form much more slow! 
than the free-fall rate in GMCs ( Zuckcrman & Evan 
[1974; Rownd & Young 1999; W ong fc Blitd 1200 

3lomon,.2004: Wu et al 2005tl . and anv success- 

ful GMC model must explain why this should be. It is 
now widely held that collapse is inhibited primarily by 
intensely supersoni c motions ( Vazaucz-Scmadeni et aL, 
UpOS"; "Mac Low & Klcsscn 2004) , rat h er than magnetic 
fields LMouschovias fc Soitzer. il976t iShu et all Il987t 
lMcKeelTl989() . (Star formation is strongly suppressed 
in low-extinction regions of molecular clouds, beyond 
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what one would expect if the star formation rate were 
simply following the column density to some power, 
sugges ting that magnetic fields may p lay a secondary 
role - lOnishi et aLlJTgM 1799^ 12001 ITohnstone et a] J 
120041 - though see iHatchell et aLl 120051 who find that 
low column densities reduce but do not completely 
prevent star formation.) Simulations and analytic 
theory indicate that the observed level of turbulence in 
GMCs is s ufficient to produce t he observed rate of star 
formation ijKrumholz fc McKeel [20051. 

However, undriven supersonic turbulence decays via 
radiation from isothermal shocks with an e-folding time 
of roughly o ne clo ud crossing time (Mac Low ct al. 1993 
Stone et alJ 119981 iMac Low! 11999: .Padoan fc Nordlundl 
19991) . soundriven turbulence alone is not suflicient to 
prevent global collapse. Instead, GMCs must either be 
destroyed before their turbulent motions decay, or the 
turbulence must be continually driven. The mode of de- 
struction is intimately related to the clouds' dynamical 
state. Unless a cloud is destroyed all at once, any in- 
ternal agent of destruction is also an internal source for 
turbulent ener gy - and one st rong enough to balance tur- 
bulent decay l)Matzneill2002|) . Destruction from within 
therefore favors models t hat achieve en ergetic and dy- 
namical equilibrium (e.g. lMcKeelll989D . if only briefly. 
The alternative - that clouds are disrupted entirely by 
outside agents (e.g. Bonnell ct al. 200(1) - requires most 
of the cloud mass to remain gravitationally unbound, as 
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bound regions rapidly collapse to higher density and be- 
come impervious to external forces. This is very difficult 
to reconcile with observational estimates of GMCs life- 
times and ratios of kinetic to potential energy. We cri- 
tique the hypothesis of unbound clouds in greater detail 
in §0 

In the following sections, we investigate the properties 
of molecular clouds both stirred and destroyed by HII 
regions within the cloud volume. The models we present 
here are improved in several ways relative to prior work: 

- Rather than enforcing strict mechanical or ener- 
getic equilibrium, we solve for the time evolution 
of a cloud's radius and turbulent velocity dispersion 
according to the time-dependent virial and energy 
equatio ns. In an early discussion of this problem, 
iMcKeel fl989f) allowed for time dependence in the 
energy equatio n, b ut assumed viria l equilibrium. 
iMatzned l|l999(l and lMatzner &: McK cc (1999) fol- 
lowed the evolution of GMCs using time dependent 
virial and energy equations, but neglected several 
terms in these equations that we model here. While 
our approach is still not a full numerical solution of 
the equations of gravity, radiation, and magneto- 
hydrodynamics, this approach enables us to study 
cloud dynamics without ignoring the effects of feed- 
back on CMC evolution, as m ost numerical sim u- 
lations to date have done fe.g. iClark et alJ|2005ft . 

- We account self-consistently for the recoil of cloud 
matter from the sites of mass ejection. In addition 
to driving turbulence, inward recoil confines the 
cloud. Recoil confinement is equivalent to an addi- 
tional (and variable) external pressure, which be- 
comes dynamically significant when cloud destruc- 
tion is rapid. Corresponding terms appear in the 
virial (§ EH Appendix Ej) and energy f H2.2l Ap- 
pendix^ equations. 

- Our model for HII regions (' ii3.2|l accounts for the 
scale-dependence of density and velocity structures 
within GMCs. Although this does not fully ac- 
count for the three-dimensional structure of the 
turbulent cloud medium, it is a significant improve- 
ment over the uniform cloud model employed by 
IWhitworth (1979). WiUiams k McKec (1997.) . and 
iMatzneri (120021) . 

- We apply the lKrumholz fc McKed ()2005|) prescrip- 
tion for the star formation rate within turbulent 
clouds. This formula accurately predicts the star 
formation rate on a variety of scales, from starburst 
galaxies to the dense precursors of individual sta r 
clusters (|Tan et al.ll200l iKrumholz fc"TMll200^ . 
We use it to govern the birth rate of ionizing asso- 
ciations (§E1). 

- Our dynamical simulations (§ track the forma- 
tion and evolution of many individual HII regions. 
This approach accounts for the finite lifetime of ion- 
izing stars and the time delay associated with the 
deceleration of shells driven by HII regions, neither 
of which is negligible compared to CMC dynamical 
times. 



We analyze the results of our models in § |H| discuss 
the implications of our findings in §[3 and summarize our 
conclusions in §|S1 In a future work (Matzner, Krumholz, 
& McKee, 2006, in preparation, hereafter Paper II) we 
apply this model to the problem of CMC formation and 
evolution in the galactic environment, including the ef- 
fects of spiral density waves. 

We do make several approximations in our work fS I7.4|l . 
We assume that the clouds are spherical and that they 
expand or contract homologously. We assume that the 
clouds are sufficiently massive that the energy injection 
is dominated by HII regions, not protostella r outflows, 
which based on the models of lMatzneil ()2002D should be 
true for clouds of mass ~ 10^ Mq or more. Such clouds 
contain most of the molecular mass in the Galaxy. We ne- 
glect the possibility that the column density of the cloud 
must exceed a threshold in order for stars to form (e.g 
IMcKeell989l) . We neglect energy injection by HII regions 
after they reach pressure equilibrium with the surround- 
ing medium. Finally, we neglect possible time dependent 
effects due to the ambient medium of the CMC: no mass 
can be added to the cloud, and the ambient pressure re- 
mains constant. Despite these approximations, the mod- 
els we present in this Paper illustrate the degree to which 
CMC properties can be understood in terms of internal 
dynamics. 

Obviously, real GMCs are not homologously expanding 
or contracting spheres with smooth density distributions, 
so the approximations we make to render the problem 
analytically tractable are quite limiting. For this rea- 
son one might be tempted give up on analytic treatment 
altogether and simply attempt to solve the problem nu- 
merically. Unfortunately, full numerical simulation of a 
giant molecular cloud, including the formation of multi- 
ple star clusters and the effects of their feedback, is not 
feasible with current codes and computers. Instead, one 
is forced to make numerous approximations regardless 
of whether one takes a numerical or analytic approach. 
Many numerical simulations of CMC evolution si mply 
ignore feedback altogether (e.g. iClark et al.l 12005V in- 
clude it only from a single source and/or focus on size 
scales much smaller than an entire GMC (^e.g. lPale et all 
l200,^lLrfc Nakamural200(t) . or focus on the galactic scale 
and lack the resolution to say anything about in divid- 
ual GMCs fe.g. ISlvz et al]l2005t iT^ker fc Brvanll200l . 
For this reason, an analytic approach that allows us to 
include feedback provides a valuable complement to nu- 
merical results, and points out areas for future simula- 
tion in which new effects might be discovered by a more 
thorough treatment of the physics. 

2. EVOLUTION EQUATIONS 

We are interested in the global evolution and energy 
balance in GMCs, so we construct a simple model in 
which we neglect details of cloud structure. We consider 
a cloud with density profile p — pe{r / Rc\)~^'' , where Rc\ 
is the radius of the cloud edge, and pe is the edge density. 
As we discuss below, we take kp = 1 as typical of GMCs. 
We approximate that the cloud evolves homologously, so 
that kp is constant. However, the cloud can expand or 
contract and can lose mass (via evaporation of gas by 
HII regions), so p^ and i?ci both vary in time. The cloud 
is embedded in an ambient medium of negligible density 
and constant pressure Pamb- We model evaporation of 
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gas from the cloud as a wind, into which cloud material is 
injected at a rate p (which we take to be negative). Gas 
that is injected into the wind travels radially outward 
with "kick" velocity relative to the radial velocity of 
the cloud at that radius. Homology requires that the 
mass loss rate follow the existing density profile, so 

. Mci 



P - 



-P, 



(1) 



where Md is the total mass of the cloud. We assume that 
the wind is low density and escapes from the vicinity of 
the cloud quickly, so that we can neglect its gravitational 
interaction with the cloud. As we discuss in more detail 
in § 13.21 this is a reasonable model for mass evaporating 
from an HII region. 

We neglect the possibility that turbulent motions 
within the cloud are likely to lead to a significant loss of 
mass. This seems justified, since in a GMC with a virial 
ratio of unity, roughly the value for observed GMCs, the 
3D turbulent velocity dispersion is smaller than the es- 
cape speed from the GMC surface by a factor of 2. 
Since the distribution of velocities in a supersonically 
turbulent medium cuts of f exponentially above the tur- 
bulent velocity dispersion l)Krumholz et al . 2006a), there 
is a negligible amount of mass moving rapidly enough to 
escape. We also neglect the possibility that a GMC might 
gain mass during its evolution, due to continuing infall. 
This assertion is more problematic, and we discuss it in 
more detail in S 17.41 

Within the limitations of these assumptions, we derive 
the Eulerian Virial Theorem (EVT) and equation of en- 
ergy conservation in Appendices ^ and We then use 
these to construct evolution equations for the the cloud. 

2.1. Equation of Motion 

We derive the equation of motion from the EVT for an 
evaporating homologously-moving cloud. 



2aiMaRc\Rc\ 
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(2) 

The proof of this theorem is in Appendix ^ In this 
equation I^i is the cloud moment of inertia, T is the total 
kinetic and thermal energy, Tq is the energy associated 
with the confining external pressure, B and W are the 
magnetic and gravitational potential energies, and the 
surface integral represents the rate of change of the flux 
of inertia across the surface iSvir that bounds the volume 
to which we apply the virial theorem. These are all terms 
that appear in the EVT for a cloud without a wind, and 
their formal definitions are given in the Appendix. The 
three additional terms represent the second derivative of 
cloud inertia caused by mass loss through the wind (the 
first two extra terms) and the rate at which recoil from 
the process of launching the wind injects momentum into 
the cloud (the final additional term). 

We now evaluate each of these terms in the context 
of our model. The moment of inertia is Id = oiMdi?^), 
where ai = (3 — kp)/{5 — kp), so its second derivative is 

-ic\ ~ a\Mc\Rl\ + a\Mc\Rc\Rc\ 
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Next consider the kinetic term T, which we evaluate 
by decomposing the velocity into large-scale homologous 
and fluctuating turbulent components (equation lASp . 
This gives 

r = ^AfelC^i -I- iaiMeii?2j + Ttu^b + 27rPamb(i?vir " ^cl)> 

(4) 

where Cd is the sound speed in the cloud (assumed con- 
stant), T^urb is the term for turbulent motions, and the 
last term comes from the constant ambient pressure Pamb 
outside the cloud. We return to T^urb below. Note that 
our assumption of homologous motion implicitly neglects 
the possibility of significant rotational motions. Ob- 
served GMCs have negligible kinetic energies in overall 
rotation compared to turbulent motions or gravitational 
potential energy. 

We use the same strategy for the gravitational and 
magnetic terms as for the kinetic term, dividing them 
into steady and fluctuating parts. The non-turbulent 
gravitational part is 
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where 
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(6) 



In principle we should include a component of poten- 
tial energy due to stars, but all observed molecular 
clouds have gas masses that greatly exceed their stel- 
lar masses. For this reason, we may neglect the stellar 
mass. For the non- t urbul ent magnetic component, we 
follow IMcKee et all Ij 19931) . Let $ be the total mag- 
netic flux threading the cloud, so that the mean field 
within the cloud \s B = $/(7ri?^j). Fr om th e the form 
of the magnetic energy term feauation lAlSp . it is clear 
that the non-turbulent component of this term scales as 

•Bnon-turb « B R^^. Wc therefore define the constant b 
such that 

$2 



B 



turb 
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The exact value of b depends on the topology of the 
magnetic field and on the background field Bq, but it 
is generally of order unity, with b = 0.3 as a typical value 
iMcKee et al.ll993D . We now define the magnetic critical 
mass Af$ by 

/ 56 \ $2 
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so that we have 



_3 GMl 
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We define the magnetic support parameter by 

A//$ 



'7b 



Afd 
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With this definition, we can combine the non-turbulent 
magnetic and gravitational terms to find 

3 GAI^ 

W„o„-turb + i3„on-turb = --a(f - vl) ■ (H) 

5 i?d 
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Now consider the turbulent components. First, we 
can neglect the turbulent component of the gravita- 
tional term because since most sub-regions of a molec- 
ular cloud are not self-gravitating (Krumholz & McKe^ 
12005*) ■ and therefore have negligible potential energy in 
comparison to their kinetic or magnetic energies. We 
can therefore set W — Wnon-turb- For the magnetic 
and kinetic turbulent components, iMcKee fc Zweibel 
(p992) ar gue for equ i partit ion of kinetic and magnetic 
energy. iStone et a l. ( 1998) find in simulations of low 
plasma-/? turbulence that magnetic energy is sl ightly sub- 
equipa rtition, ;Bturb ~ 0.6 7^urb, which M cKee fc TanI 
argue can be understood as the kinetic and mag- 
netic energies reaching equipartition for motions trans- 
verse to the field, but not for motions along the field. 
W e adopt the ratio o f magnetic to kinetic energy found 
bv IStone et al.l l)1998j) . so the combined turbulent kinetic 
and magnetic energies in the virial theorem are 
2rturb + Sturb - 2.6Tt„b - "2 

where ad = ((w^ - Uci,^)^)^^ 
mass-weighted turbulent velocity dispersion of the gas 
in the cloud. 

Finally, we come to the surface terms, Tq and 
(d/dt) Jg [pvr^) ■ dS. We can make the latter term zero 
by choosing our virial surface to be well outside the cloud 
so that the density of cloud material on the surface is neg- 
ligible. However, the pressure outside the cloud Pamb is 
non-zero, so 



3.9Mdati, (12) 
is the one-dimensional 
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Substituting into the EVT ((2Jl, we arrive at an equation 
of motion for the cloud: 

aii?ci = 3.9^ + 3|L-|(l-^i)a^ 



amb 



cl 



VI 
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cl KpJ Mcl 

This equation is intuitively easy to understand. The left- 
hand side represents the acceleration of the cloud edge, 
which is equated with the force per unit mass due to 
internal turbulence and pressure (the first two terms), 
gravity and magnetic fields (the third term), external 
pressure (the fourth term) , and recoil from the evaporat- 
ing gas (the final term). 

For convenience we wish to non-dimensionalize this 
equation. Let Mci_o, i?ci-o, and Uci-o be the ini- 
tial mass, radius, and velocity dispersion of the cloud. 
We define the dimensionless variables M — Md/Md-o, 
R = Rci/Rci-Q, cr = CTci/o'ci-o, and r = i/^ci-o, where 
^cr-o ~ Pci-o/o'ci-o is the crossing time of the initial 
cloud. In these variables H14|l becomes 
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where the primes indicate differentiation with respect to 
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(16) 
(17) 
(18) 



and we have defined the initial Mach number 

O'cl-O 



Ccl 
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and non-thermal virial parameter l)Bertoldi fc McKe^ 
IT99I 



Ckvir-O — 
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(20) 



If the cloud's initial state is in equilibrium and it is not 
losing mass, so R'{0) = M'(0) = 0, then the ambient 
pressure Pamb must be such that 



3.9 + 37^0^ 

VP = VG- 

ai 



(21) 



Note that this will generally implies an ambient pressure 
higher than the mean in the ISM. We consider it realistic, 
however, because GMCs form in relatively overpressured 
regions, and because we have not included the weight 
of an atomic layer overlying the cloud. Once mass loss 
commences, there will be an additional recoil pressure. 

2.2. Equation of Energy Evolution 

To derive the evolution equation for the cloud energy, 
we begin with the general energy conservation equation 
for an evaporating homologous cloud, which we derive in 
Appendix IbI and for convenience repeat here: 

f = ^ [£ + (1 - ^i)W] - 4.Pa.„bP?,i?cl 



3-fc, 
A-k,, 



MclPclfoj + Gc\ - Cc\. 



(22) 



Here £ is the total cloud energy, and Gci and Cc\ are the 
rates of radiative energy gain and loss integrated over the 
entire cloud. This equation is easy to understand intu- 
itively. The term {Mci/AIci)£ is simply the mass loss rate 
times the energy per unit mass in the cloud. The term 
(Mci/Mci)(l — ?7b)W is the rate at which mass loss re- 
duces the cloud energy via reduction of the gravitational 
and magnetic fields, both of which are proportional to the 
mass. The next two terms represent the rate at which 
the external pressure and the recoil force from launch- 
ing the wind do work on the cloud. Finally, the last two 
terms are simply the rate of radiative gains and losses. 

Using the same arguments as in § 12.11 we may write 
the total energy in the cloud as 



aiMdRlx + 2AMd(jly 



Mdcly 



7«(1-%)^T^- 

5 Rc\ 

(23) 

Note that the factor of 2.4 in the Mcia^i term comes from 
taking T^urb + '^tmb ~ l-6'7turb, and the 3/2 in front of 
the Mc\c^\ term comes from the assumption that the ra- 
tio of specific heats for the cloud is 7ci = 5/3. One 
might expect 7/5 instead, since the cloud is molecular 
and therefore diatomic. However, the lowest rotational 
or vibrational excitations of H2 have excitation temper- 
atures of several hundred K. Since the cloud is far colder 
than this, molecules never have enough energy to access 
their rotational and vibrational degrees of freedom, and 
the gas acts as if it were monatomic. The time derivative 
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of this is 

d£ 1 ■ ■ n • •• • n 

— = -aiMciRli + aiMciRc\Rc\ + 2.4 M^al 



4.8 AfciCTciCTci^MciCci - -a(l - 77I 



3.-. 
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Rc 



(24) 



Substituting £ and d£/dt into the energy equation (|22|l . 
we find 

aiMcli?cl-Rcl + 4.8MclCrclCTcl 

-47rPambi?cl4l + ( ) Mcli^cll'ei 



(25) 



We may regard this as an evolution equation for a^i, 
which makes intuitive sense: the overall expansion and 
contraction of the cloud is dictated by the equation of 
motion, and the thermal energy per unit mass is fixed, 
so the turbulence acts as the energy resevoir, increasing 
or decreasing as the cloud gains or loses energy. Re- 
arranging to solve for &ci and non-dimensionalizing as 
we have done with the equation of motion gives 



ai 



R'R" 



AIR' 
R^' 



-rip- 



R'R' 
Ma 



-VR- 



M'R' g~c 



Ma aiMa ' 
. (26) 

where G and C are the dimensionless rates of radiative 
energy gain and loss, defined by 



g-c 



Rci-0 



(^c 



(27) 



3. ENERGY SOURCES AND SINKS 

In this section we evaluate the rates of radiative en- 
ergy gain g and loss C, and the characteristic launch 
speed of evaporating gas . Together with the equation 
of motion H15(l and the energy equation H26|l . and the 
star formation rate (§0}, this will completely specify the 
evolution of our model clouds. 

3.1. Decay of Turbulence via Isothermal Shocks 

GMCs are approximately isothermal because their ra- 
diative time scales are much shorter than their me- 
chanical time scales. As a result, supersonic motions 
within the cloud generate radiative shocks that remove 
energy from the cloud. The problem of the decay of 
turbulent motions by supersonic isothermal shocks in 
both hydrodynamic and magnetohydrodynamic media 
has been studied ex tensively by nume rical simulation 
(e.g. 'Mac Low et al. 1998; Stone ct al. "1998"; 'Mac Low| 
0999; Padoan fc Nordlund 1999). Stone et al. (199§ 
finds that the dissipation time scale of the turbulent en- 
ergy is tdis = E/E = 0.83Ain/wrms = 0.48Ain/cr, where 
Ain is the length scale on which the energy is injected. 

In reality, HII regions coming from associations of var- 
ious sizes, winds, and gravitational contraction of the 
cloud will all contribute to turbulent motions and inject 



energy on different length scales. However, we can es- 
timate an effective length scale through a combination 
of observational and theoretical considerations. Obser- 
vationally, turbulence in nearby GMCs appears to be 
driven on s cales comparable to the cloud scale or larger 
(iBasu fc Murah 2QQI: ..Hever fc Br unt 20041 and theo- 
retical estimates of the effective energy injection scale of 
HII regions suggest that this is also near the cloud scale 
(Matzner 2002). The longest wavelength mode that a 
cloud of radius i?ci can support is ~ 4_Rci, where the 
factor of four arises because the largest turbulent mode 
corresponds to overall expansion or contraction of the 
cloud, in which diametrically opposed points are mov- 
ing in opposite directions. Motion in opposite directions 
corresponds to the points being half a wavelength apart, 
lying a total wavelength of twice the cloud diameter 
Matzneil I2OO2') . Thus, we take the effective injection 
scale to be Ain = 4(/>ini?ci, where (pin < 1, and we take 
(pin = 1 as a fiducial value based on observations and the- 
ory. The (dimensionless) rate at which energy is radiated 
away due to decaying turbulence is 



?7v Ma^ 

(bin R 



(28) 



where Stone et al.'s simulations give rj^ — 1.2. Note that, 
in deriving this factor, we have used our result that the 
energy in the turbulence, including magnetic and kinetic 
contributions, is 2.4Mcicr^j. 

It is also worth noting the poss ibility that the me asured 
energy loss rates are too high. iGho fc Lazari an (200^^ 
argue that Alfven waves cascade and decay anisotrop- 
ically, and that this anisotropy can reduce the decay 
rate. However, Cho & Lazarian argue that simulations 
to date have not captured this effect bec ause they use an 
isotrop ic driving field that is unrealistic. iSugimoto et alJ 
lj2004[l simulate filamentary molecular clouds, and find 
that Alfven waves of certain polarizations decay more 
slowly than simulations have found. Even if neither of 
these effects apply, there are differences in the rate of de- 
cay in different simulations depending on how the tu rbu- 
lence is driven. The simulations of Mac Low! l)1999D give 
slightly lower dissipation rates, probably because in those 
simulations the turbulence is forced with a driving field 
that is constant in time, while in those of Stone ct aQ 
((T998() the driving field is determined randomly at each 
time step. While we feel that the Stone et al. approach 
is somewhat more realistic, this is by no means certain. 

3.2. HII Regions 

HII regions are the do minant source o f energy injection 
into GMCs from within l)Matzneil2002() . We consider the 
effects of the HII region from a single association here, 
and extend our results to a population of associations to 
§11 

3.2.1. Evolution of Individual HII Regions 

We first describe the evolution of an HII region em- 
bedded in a mo lecular cloud, modifying the analysis by 
iMatzned l)2002|) by allowing the mean ambient density 
and velocity dispersion to vary with radius r away from 
the formation site of an association, as p{r) oc r~^'' and 
a[r) (X r"*^", respectively. iMatzneJ l)2002j) considered 
only the homogeneous case kp = ka = 0. Note that the 
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local turbulent virial parameter a{r) = 5ra'^{r)/[GM{r)] 



scales as 



Since a{r) is roughly unity on the 



scale of the cloud and on that of the newborn associa- 
tion (else it would not have formed), it is reasonable to 
assume ko = 2ka- + 2. The observed line widt h-size rela- 



tion and density-size relations l|Solomon et ai.iiitjSYii im- 
ply kcr = —1/2 and fcp = 1, and we adopt these values be- 
low. These parameters correspond to a cloud with a neg- 
ligible internal pressure gradient. Note that it has been 
suggested that the observational result kp = 1 (which 
is equivalent to GMCs having constant column densities 
within a galaxy) is simply an o bservational artifact (e.g. 
iVazauez-Semadeni et al.lll997|) . However, many of the 
proposed mechanisms to explain how this artifact could 
be created do not apply to extragalactic observations, 
and more recent observations show that GM Cs in other 
galax ies also show constant column densities (jBlitz et al.l 
120051 and references therein). We discuss this point in 
more detail in S 16 .21 and S I7.3| and also refer readers to t he 
discussion of this point in Krumholz fc McKed l)200fiD . 

The density variation affects the expansion phase of 
the HII region, and the variation in velocity disper- 
sion affects how HII regions merge with the background 
turbulence. The evolution of an HII region in a tur- 
bulent CMC is a substantial problem that must ul- 
timately be solved via simulation. However, to date 
only preliminary attempts to solve the problem have 
been made (e.s.. [Dale et al.l 120051 iMellema et all 120051 
iMac Low et al.ll2006HKrumholz et alJl2006b|) . so we are 
forced to rely an analytic approximations. 

First, consider the expansion phase of an HII re- 
gion. Assume that it has expanded well beyond its 
initial Stromgen radius. The mean ambient density is 
p(r) oc r~'"'p within a distance r from the association, 
the ionized gas temperature is Tn — 7000 K, the (con- 
stant) ionizing luminosity is S* = 10"*^ 6*49 s~^, and the 
recombination coefficient is a*^^^ in the on-the-spot ap- 
proximation. If the HII region is blister-type, ionized 
gas will rocket away from the ionization front at a ve- 
locity tigj = 2cii, where cn = 9.74 km s~^ is the sound 
speed in the ionized gas. This is the characteristic launch 
velocity of our cl oud's esca ping wind. For the evolution 
of the HII region, 'Matznei*s dynamical equation [his eq. 
(15)] admits the self-similar solution 



rLpirsid = (2, 1) X - 



3 (7 2fcp)^ ^2 



4 9-2fco 



-PllCiit 



(29) 



for (blister, embedded) regions, respectively. Here rgh is 
the radius of the shell at time t, and pu is the (uniform) 

3/2 

density inside the HII region, which must vary as r^j^ 
to ensure that ionizations balance recombinations. This 
equation applies when the expansion velocity is much 
larger than the turbulent velocity dispersion; we discuss 
the effects of turbulence below. This yields 

rlCpir^^) = {2,l)x2.2^g^^ksT,/ 



2(9-2/cp) 



\47ra(2) 



(30) 

where ks is Boltzmann's constant. The leading factor of 
2.2 accounts for the particle abundances assuming helium 
is singly ionized. We assume blister-style regions here- 
after, because GMCs are porous and therefore even small 
HII regions are likely to be able to punch through to low 



density channels and vent l)Matzneill2002|) . For our fidu- 
cial case fcp = 1, the mean column S = 4p(r)r/3 is a con- 
stant. Adopting a'^^ = 3.46 x 10~^^ cm"^ s~^, and adjust- 
ing the effective ionizing luminosity d ownward by a factor 
1.37 to account for dust absorption ijWilliams fc McKed 
[1^3): we find 



r,h = 15.451/' 



t 



3.8 Myr 



4/5 



0.03 g cm 



2/5 



pc. 
(.31) 

Here we have normalized E to a typical value for Milky 
Way molecular clouds. We take the typic al ionizing life- 
time to be 3.8 Myr because we adopt the Krouoa (2001]) 
initial mass function (IMF) and the st ellar ionizing lu- 
minos ities and lifetimes tabulated by iParravano et alJ 
l)2003(l . However, the lifetime can fluctuate from this 
for small clusters, so we define as the main-sequence 
ionizing lifetime of stars in a given cluster. (With these 
choices, a young association large enough to fully sample 
the IMF emits 3.4 x 10^^ ionizing photons per second 
per solar mass, and has one star above 8M0, hence one 
supernova, per 131 solar masses.) The momentum of the 
HII region during this driving phase is 



fth = 5.1 X 1043^fg/5 



t 



0.03 g cm 



3.8 Myr 

-2\ 1/5 



7/5 



g cm s 



(32) 



Assuming the gas can escape freely, the mass evaporation 
rate from the cloud is simply 



-]3sh/(2cil). 



(33) 



Note that equations H31|l and (|32|) are quite insensitive 
to the actual escape of ionized gas (Matzner 2002), but 
Mci depends on the existence of an escape route. 

After the driving stars burn out at t = tms (and if it 
has not yet merged with the cloud turbulence, ? I3.2.2|l . 
the HII region will continue to expand in a momentum- 

— (3 — ) _2 

conserving snowplow, in which fgh r^^ '' r^j^ so 
that 



rsh 



1/3 



(34) 



where r^s and r^g are the radius and velocity of the 
expanding shell at the time when the driving stars burn 
out. No mass is evaporated in this phase; we ignore the 
possibility of dynamical mass ejection by the snowplow 
unless the snowplow contains enough kinetic energy to 
unbind the cloud as a whole. 

We now comment on two potential criticisms of our 
HII regions. First, real molecular clouds are very in- 
homogeneous, being filled with clumpy and filamentary 
structure. To what degree should th is affect the re s ults o f 
this section? Very little, we argue. IMcKee et al.1 Ipsl 
consider HII regions in a medium composed entirely of 
dense clumps. Adopting a clump mass distribution very 
similar to that observed within molecular clouds, they 
find (their eq. 5) that the photoevaporative rocket effect 
clears all but a couple clumps from withi n raVi(t). Addi- 
tional homogenization, not considered bv iMcKee et all 
should come from the ram pressure stripping of overdense 
clumps and filaments as they are struck by gas already 
set in motion. For these reasons we expect the gross 
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properties of HII regions to be rather ins ensitive to local 
inhoin ogen eities. Recent s i mulat ions by iMellema et alJ 
()2005|) an dlMac Low et all l)200()|) support this expecta- 
tion, and iMellenm^r'alJ find that the radius of an HII 
region in a turbulent medium as a function of time is 
very similar to what one wou ld find in a uniform medium. 
(The simulation bv lDale eTa l. (2005) gives a somewhat 
different result, but its different initial conditions and 
brief d uration do n ot allow a fair comparison. In par- 
ticular, ID ale et al.l compare their HII region expansion 
rate to the analytic solution for a uniform medium, de- 
spite the fact that they are simulating a strongly cen- 
trally condensed gas cloud.) For these reasons we expect 
our results to be robust against clumping of the gas on 
scales smaller than rsh(i)- Given that our adoption of 
blister-type flow reflects the existence of inhomogeneities 
on scales larger than rsh{t), we expect the results in this 
section to be robust against cloud inhomogeneities in 
general. 

A second potential criticism concerns our assumption 
that each HII region expands as if the gas density drops 
off like r~'''' away from it. Formally this can only be 
consistent with our approximations of a spherical, p oc 
r^'^'p cloud if every HII region is born at the cloud center. 
We consider this to be merely a technicality. While our 
global cloud model refers only to a mean density profile 
and a turbulent energy, our model for the local density 
enhancement reflects the fact that clusters are born at 
the peaks of the turbulent density profile - which will 
not always occur at the cloud center. Indeed, we model 
mass loss and recoil pressure as if HII regions peppered 
the entire cloud and commonly vent out its side, and we 
ignore the interaction between regions. 

Although our local kp = 1 density is an idealization, 
we consider it an improvement over previous works which 
assumed homogeneous cloud gas. Moreover, the most 
luminous HII regions, those which dominate the energy 
injection, are in reality likely to be born near cloud cen- 
ters, and once they have expanded to a radius that is 
significantly larger than the distance from their center 
to the peak of the density distribution in the GMC, the 
approximation that they are expanding down a fcp = 1 
density gradient will be quite good. 

3.2.2. Merging of HII regions 

The shell will continue to expand until it decelerates 
to the point where the expansion velocity is compara- 
ble to the current velocity dispersion in the GMC, at 
which point the shell will break up, merge with the gen- 
eral turbulent motions, and cease to be a coherent en- 
tity. The HII region experiences a velocity dispersion 
cr(r) = a{r / Rci)^^" after expanding to radius r, and we 
approximate that the merger occurs when (T(rsh) = Tsh- 
If the cloud properties were to remain fixed during the 
expansion, the merger radius of an HII region would be 



l/(3-fe„-fep)^2/5 



•cl, 



(35) 



2psh 

applying fcp ^ 1 and ka 
from our idealization of the cloud as a sphere and the 
HII region as a hemisphere. Equation (|35|l holds whether 
merging occurs in the driving phase (psh ^ Psh{t)) or the 
snowplow phase {psh Psh(ims))- However, since the 



-1/2. The factor of 2 arises 



cloud velocity dispersion and radius can change during 
the expansion of an HII region, equation H35(l holds only 
approximately. 

Note that gravity has little effect on HII regions ex- 
cept during merging, because their shells travel faster 
than a{r) up to that point. Indeed, the crossing time 
r jair) is 2.0a(r)^^/^ times the free-fall time on scale r. 
Since a(r) > 1, gravity only marginally affects shell mo- 
tions. For this reason, we are justified in continuing to 
use our similarity solution (which neglects gravity) up to 
the point at which a shell merges. 

3.2.3. Energy injection 

At any given time, the cloud energy terms T and Ai 
include both the turbulent energy from merged HII re- 
gions, and the kinetic energy of those that are still ex- 
panding. During its expansion, the kinetic energy of a 
single HII region is 



^ 1 

^1 = ^^sh^'sh 



(36) 



the two terms reflecting bulk and thermal energy, re- 
spectively, if Pii is the internal pressure. The bulk 
kinetic energy at the time of merging tm is therefore 
Ti{trn) = Pshcr{rm)/2. 

Consider flrst the energy injected into turbulence at 
merging. We assume a merging shell stores the same 
fraction of energy in non-kinetic form as does the back- 
ground turbulence, so we consider the total energy con- 
tribution to be 1.67i(tm). Rather than track individ- 
ual regions' contributions in detail, we lump them to- 
gether in the energy budget in a way that reproduces the 
correct average energy in turbulence. Since the cloud's 
energy dissipation time is tdis = 2.O(1.20in/77v)i?ci/o'ci, 
whereas the dissipation of a single contribution occurs 
over tdis,i = 2.O(1.20i„/77v)r„/cr(r™), we add 

1/2 



AE< 



turb.l 



X l.eTi ^ ?7E (^-^ j xl.eTi, 

(37) 

rather than 1.67i, to the turbulent energy. The factor 77E 
parameterizes our ignorance of the exact efficiency with 
which HII regions drive motions in molecular clouds. For 
a uniform medium t^e = 1, but there has been some sug- 
gestion that tik < 1 in turbulent media (e.g. tElmcgrce^ 
I2OOOI: IPale et al.1 120 05"). although other simu lations do 
not appear support this result fe.g. IMellema et 31,^200^] . 
In our models we vary t/e to see explore how sensitively 
our results depend on the efficiency with which HII re- 
gions inject energy. 

We exclude from consideration the energy that expand- 
ing shells possess prior to merging (which we estimate to 
be several times smaller than the mean turbulent energy 
contributed by these regions). We note that the cloud 
outside of a given region is unaffected by it, except grav- 
itationally, until it has merged. Thus, the rate of energy 
injection due to a single HII region is roughly 

1/2 

6{t ~ t„0, (38) 



Qc\ = 1.677e^i 

and when a shell merges with the overall cloud turbulence 
we set 
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1/2 



(39) 
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3.2.4. Cometary HII regions and cloud disruption 

Finally, note that it is possible that a shell will con- 
tain enough momentum to expand to the point where 
its radius is larger than the cloud radius. If the expan- 
sion velocity of the shell at this point is larger than the 
escape velocity of the cloud (i.e. Tsh > \/'2GMc\/ Rc\), 
then the shell will simply unbind the cloud. This crite- 
rion for cloud disruption is somew' h at diff erent fr om those 
adopted by iWilliams fc McKed 1)1997(1 and iMatzneJ 
(j^02[); we discuss the distinction and its implications in 
S I6.1I If the expansion velocity is smaller than the escape 
velocity, the s hell will deform the cloud into a cometary 
configuration l)Bertoldi fc McKe'^ll99(]|) . but what hap- 
pens at that point depends on whether or not the shell is 
still being driven by an active association. An undriven 
shell will neither gain nor lose energy once rgh > ^ci, 
and its energy will eventually be added to the cloud's 
as turbulence once the cloud falls back. We therefore 
approximate that any shell in the snowplow phase that 
reaches rgh > Rd but does not have enough momentum 
to unbind the GMC merges with the turbulence. On the 
other hand, a shell with an active source can continue to 
gain energy and evaporate mass even after reaching the 
cometary configuration, altho ugh eventually its en ergy 
will saturate. Following Williams fc McKea l(1997tl . we 
estimate that a shell can continue gaining energy after 
reaching the cometary phase up to a time t ^ 2ti^, where 
tn is the time it took the HII region to reach Rc\. Its 
radius at this time will be Tsh = 1.74i?ci. If a driven HII 
region reaches this radius and still does not have enough 
energy to unbind the cloud, then it can affect the cloud 
no further. We approximate that the HII region merges 
with the turbulence and evaporates no further mass after 
that point. 

It is important to point out that our treatment of large 
HII regions, and our criteria for cloud disruption and 
deformation into a cometary configuration, are quite ap- 
proximate. In comparison to lMatzneil l)2002fl . the criteria 
we use here reduce the shell momentum required to dis- 
rupt a cloud by roughly a factor of 2. This may affect 
our conclusions about the frequency of disruption by HII 
regions. 

3.3. Supernovae and Protostellar Outflows 

Here we discuss in more detail two sources of feedback 
that we have chosen not to include: supernovae and pro- 
tostellar outflows. We omit these because they deliver 
much less energy per unit mass of stars than HII regions. 
Here we summarize the arguments, many of which are 
given injKlatzner (2002), as to why this is the case. 

The dynamics of a supernovae exploding inside HII re- 
gions, and the amount of energy and momentum they 
add, may be studied directly either by numerical simu la- 
tion fe.g. lTenor io-Tagle et al."1985':' Yorke et al.lll989n or 
analytic calculation (Matzner 2002). Both methods yield 
similar results: supernovae generally sweep up enough 
mass inside their ionized bubbles to become radiative by 
the time they reach the outer boundary of their host HII 
regions. As a result, supernova blast waves radiate away 
much of their energy before encountering any molecular 
cloud material, and consequently increase the total mass 
removed by their HII region by at most ~ 20 — 40%, and 
more typically ^ 10%. Moreover, this applies only to 



supernova progenitors whose HII regions are enclosed by 
the cloud and do not blister. Supernovae inside blister- 
type regions will deposit most of their energy outside the 
cloud, into the low density ambient medium, and should 
therefore have an even smaller effect. The as sumption 
sometimes made in the literature fe.g. lClark et al. 200ftl 
that supernovae can halt star formation and disrupt 
molecular clouds in ^ 10 Myr does not appear to be 
supported by any numerical or analytic calculations. 

Protostellar outfiows are negligible compared to HII 
regions on the size scale of entire GMCs for two rea- 
sons. First, the momentum injected into a cloud by an 
outflow is of order M^,v^, where is the star forma- 
tion rate and Vw ~ 40 km is the rough momentum 
input per unit mass of stars formed for hydromag netic 
winds l|Matzner fc McKedl2000l iRicher et al.ll200o|) . In 
contrast, as the analysis of § 13.21 shows, the momentum 
carried by HII regions is of order AfphotCii, where Mphot is 
the rate at which ionizing photons photoevaporate cloud 
mass. While Vw is larger than cn by roughly a factor of 4, 
the low star formation efficiencies of GMCs imply that 
the rate at which ionizing photons evaporate gas must 
exceed the rate at which the gas transforms into stars 
by a factor of 10, so Mphot ^ M^,. Our results in §|H1 
confirm this conclusion. Thus, HII regions provide much 
more momentum than outflows. 

Second, outflows generally inject energy on signifi- 
cantly smaller scales than HII regions. While there are 
some spectacular examples of I IH objects parsecs in size, 
these appear to be exceptions. lOuillen et al. (2005) find 
in NGC1333 that the typical size of protostellar out- 
flow cavities is only ~ 0.1 — 0.2 pc. This is a low- 
to intermediate-mass star forming region, and cavities 
are likely to be even smaller in the denser environments 
where most Galactic star formation takes place. This 
is much smaller than the sizes comparable to the cloud 
radius reached by large HII regions (Matzner 2002). Fol- 
lowing the discussion in § 13.11 we expect the energy in- 
jected by these outflows to decay much more quickly than 
the large-scale motions created by HII regions, and there- 
fore make an even smaller contribution that comparison 
of the momenta would suggest. 

It is worth noting that both arguments show that HII 
regions dominate over outflows only apply if we are con- 
cerned with objects comparable in size to entire GMCs. 
For smaller objects like the gas clumps forming individ- 
ual clusters, the time for a single HII region to expand 
may be comparable to the evolution time of the entire 
region, and much of the energy of the HII region may be 
deposited outside the clump. As a r esult, outflows may 
well do minate for these objects (e.g. iQuillen et al.ll2005t 
iLi fc N akamura 2006). 

4. STAR FORMATION AND HII REGION CREATION 

To know how HII region feedback affects GMC evolu- 
tion, we must know the star formation rate and the clus- 
tering properties of the stars formed. To estimate the 
former, we use the turbulenc e -regu lated star formation 
model of lKrumholz fc McKe j l|2005|) . This model postu- 
lates that stars form in any sub-region of the cloud where 
the local gravitational potential energy exceeds the tur- 
bulent kinetic energy, and it gives good agreement with 
simulations, with the observed star formation rate in the 
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Milky Way llLuna et al.ll2006t), w ith the age spreads of 
rich star clusters ijTan et alJl200fiD . wi th the star forma- 
tion rate in dense molecular clumps (^Krum holz fc TanI 
I2006I) , and with the the Kenn icutt-Sch midt Law for star 
formation in galactic disks l)Kennicu tt 1998). In the 
Krumholz & McKee model, the star formation rate is 



A'U = SFRff 



Mel 

iff 



(40) 



where 



SFRff' « 0.073a~°-^^M 



(-0.32 



(41) 

iff = [37r/(32G'p)]i/2 is the free-fall time at the mean 
density p of the cloud, and Ovir and M are the virial pa- 
rameter and ID Mach number. Thus, given the mass, ra- 
dius, velocity dispersion, and gas temperature in a GMC, 
the Krumholz fc McKee model enables us to predict the 
instantaneous rate of star formation in that cloud. 

HII regions will be driven by the combined ionizing 
luminosity of all the stars in an association, and the en- 
ergy injection therefore depends on how the stars are 
clustered. Thus we must estimate the ionizing luminosity 
function of associations as well as the star formation rate. 
This distribution is unfortunately quite uncertain, be- 
cause a small cloud cannot form arbitrarily large associa- 
tions and thus the Galactic-average and individual-cloud 
luminosity functions are different. Let dFa{N'^) / dlnj\f^ 
be the fraction of OB associations with In number of stars 
between In A/"* and In A/"* -|- dlnA/"* in the entire Galaxy, 
and let di^a,Mci(A/'*)/dlnA/'* be the corresponding frac- 
tion within a cloud of mass Mci. iMcKce fc Willian3 
((1997') show that observations of the Galactic population 
of HII regions are consistent with a Galactic distribution 



1 



oc 



dlnAA* A/", 



(42) 



It will be more convenient for us to work with the mass 
of an association rather than the number of stars. Since 
we are concerned with OB associations which give rise to 
HII regions, and these have enough stars to sample the 
IMF well (A/; > 100, Zinneckcr et al. 1993), we convert 
from number to mass simply by multiplying by the mean 
stellar mass. We use the .Kroupa, (2001) IMF, which gives 
(™)iMF ~ ^■'^^ -^^Q' where ()jmf indicates an average 
over the IMF. Thus, 



dK, 



1 



cx 

dlnMn AL 



(43) 



is the Galactic distribution of association mass Ma ■ 

To go from the Galactic-average distribution to the in- 
dividual clo ud distribu t ion, w e follow Williams fc_ McKec 
l|1997() and iMatzneJ l|2002() . First, we require that 
(1) no association has a mass larger than eMd, where 
iWilUams fc McKed l|1997D estimate e = 0.1, (2) the dis- 
tribution of association masses within individual GMCs, 
when integrated over the GMC population of the Galaxy, 
gives the Galactic distribution 1)43(1 . and (3) the star 
formation rate per cloud scales with the cloud mass as 
oc M^. Using an ar gument analogous to the de riva- 
tion of equation (16) in WiUiams fc McKee (1997) and 
equation (35) in .Matzne r (2002), we derive a distribution 



din Ma °^ 1 - (eMei/M„)"-/3 



Here H{x) = (1,0) for (x > 0,a; < 0) is the Heavi- 
side step function, and the Galactic population of star- 
forming GMCs satisfies dAfd/dlnMd oc M'" with 
a 0.6 to an upper l imit of M„ = 6 x 10^ Mq 
((Williams fc McKee 19 9-j). Based on th eir star forma- 
tion model, .Krumholz fc McKed l(2005|) suggest f3 w 
0.67, and we adopt this value in our work. Note that /3 < 
1 implies a higher rate of star formation per unit mass 
in smaller clouds, which occurs because in the Galaxy 
smaller clouds have higher densities and thus shorter 
free-fall times. Equation 1(44(1 is fairly simple to under- 
stand intuitively: small clouds cannot make arbitrarily 
large associations, hence the factor H{eMci — Mg). If 
all clouds produced OB associations with the same mass 
distribution as the Galactic powerlaw distribution, up to 
the maximum mass they could produce, there would be 
a deficit of large associations because all clouds can make 
small associations, but only a small fraction of clouds can 
make large ones. To offset this effect and produce the 
Galactic distribution of OB association masses, GMCs 
must be slightly more likely to produce an association 
near their upper mass limit than a straightforward ex- 
trapolation of the Galactic association mass distribution 
would suggest, hence the factor [1 — [eMcx/ Ma)°'~^\~^ ■ 
For large associations, the ionizing luminosity is simply 



5*49 — 



(S49(m)) 



IMF 



Ma, 



(™)lMF 

and the main sequence ionizing lifetime is 



(ims)a 



(s49(m)i,ns(™))iMF 



(S49(»7i)) 



(45) 



(46) 



IMF 



where 549(771) and tmsijn) are the ionizing luminosity (in 
units of 10^^ photons s~^) and main-se quence lifetime of 
a star of mass m. We adopt the fits of iParravano et till 
i2003ll f or s(m) and tins('7i), which, together with the 
IKroupa) pOOl) IMF, give (s49(m))jMF = 7.2 x 
S'49 = 3.4x 10~^(Ma/MQ), and t^^^ = 3.8 Myr. However, 
for smaller associations the ionizing luminosity is likely 
to be dominated by the single most massive star, and this 
causes the ionizing luminosity function to flatten and the 
ionizing lifetime to vary with S49. McKee fc Williani^ 
11.997, Appendix A) give an analytic formula that ap- 
proximates the flattening, but for semi-analytic models 
we can dispense with the approximation and determine 
the luminosity function simply by drawing stars ran- 
domly from the IMF until we accumulate mass up to a 
given association mass. We can then determine the ion- 
izing luminosity of the association by summing those of 
the individual stars, and define the main sequence life- 
time as the time at which the stars providing half the 
ionizing photons disappear. 

N ote that in co ntrast to (Williams fc McKed l(1997fl 
an d [MatzneilpOOl . we do not impose an absolute upper 
limit of 5*49 < 490 on the ionizing luminosity of associ- 
ations. This has no effect at all on any but the most 
massive clouds, since the requirement that the associa- 
tion mass be at most 10% of the cloud mass imposes a 
limit on 6*49 that is lower than this. For the most massive 
clouds we model in § 16.11 the largest associations formed 
reach 6*49 ~ 1700, but the extremely weak dependence 
of HII region properties on S'49 (rgh oc S\q^) means that 
increasing the maximum S'49 by a factor of a few has very 
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little effect. Furthermore, due to the relative improbabil- 
ity of forming an association so close to the upper mass 
limit, even for our highest mass models the majority of 
clouds do not form associations with ^49 > 490. Thus, 
we do not expect the presence or absence of an upper 
limit on association ionizing luminosities to affect our 
results in any significant way. 

5. SEMI-ANALYTIC MODELS 
5.1. Methodology 

We now have in place all the necessary theoretical ap- 
paratus to set up our semi-analytic models. In essence, 
our model describes the evolution of a GMC using a pair 
of coupled non-linear ODEs, with added damping and 
driving terms. We integrate these equations forward in 
time in a three-step process. First, we use the current 
configuration of the cloud to compute the rate of tur- 
bulent decay (equation I28|) , the rate of star formation 
(equation I40|l . and the rate of evaporation by Hll re- 
gions (equation 1231 ■ We update Rd, Rd, crd, and Md 
using these values (equations El ^^nd . We chose our 
update time step so that no quantity changes by more 
than 0.1% per advance. Second, we update the state of 
HII regions as described in § 13.21 We compute new val- 
ues of the radius and expansion rate for each HII region, 
removing those whose expansion rates are low enough or 
radii are large enough for them to merge. At merging, 
we add their energy to turbulent motions in the cloud. 

Third, we create new HII regions. To do this we gen- 
erate a random mass for the next association, chosen 
from the distribution (|44|) . We track the amount of mass 
transformed into stars since the last assocation was fully 
formed, and when half the next association mass has been 
accumulated in new stars, we add an HII region for that 
association. To determine the properties of the HII re- 
gion, we generate stellar masses from the Kroupa (2001D 
IMF, assign an ionizing luminosity and main sequence 
lifetime to each star, and use these to compute the to- 
tal ionizing luminosity and lifetime (defined as the time 
when the stars responsible for half the ionizing photons 
burn out) for the association. We continue to put newly 
formed stars into the new association until the full mass 
of the association has been accumulated, at which point 
we re-set the tally of accumulated stellar mass and ran- 
domly generate a new mass for the next association. 

We terminate the evolution when one of three condi- 
tions has been satisfied: (1) an HII region unbinds the 
cloud (i.e. rsh > R and fsh > y/2GM/R); (2) the cloud 
surface density has dropped to the point where its visual 
extinction Ay < 1.4, the minimum required for CO to re- 
main molecular ( van Dishoeck & Bla ck 1988), assuming 
the standard Milky Way interstellar UV field and dust to 
gas ratio, whereby 1 g cm~^ corresponds to Ay — 214; 
or (3) the time step is less than 10~* times the current 
evolution time, which occurs if the radius is approaching 
zero. We term these possibilities disruption, dissociation, 
and collapse. An important caveat that applies to all 
these outcomes is their dependence on our assumption 
of spherical symmetry. Even if an HII region delivers 
an impulse capable of unbinding a cloud, the cloud may 
actually be displaced as a whole, or it may break into 
multiple pieces, each of which is internally bound and 
capable of continuing to form stars. In the case of dis- 



TABLE 1 
Fiducial parameters. 



Parameter 


Value 


"vir-O 


1.1 


Cs 


0.19 km s-i 


Vb 


0.5 


m 


1.0 




1.2 


kp 


1.0 


4>in 


1.0 



sociation, a CMC's mean column density may be so low 
that much of its mass is turned atomic by the interstellar 
UV field, but overdense clumps within it may survive and 
continue star formation. Finally, the collapse case would 
probably result in a cloud fragmenting into smaller pieces 
rather than undergoing monolithic collapse as occurs in 
our one-dimensional models. 

5.2. Fiducial Initial Conditions 

Here we describe a basic set of initial conditions for 
our runs. In § |^ we discuss how varying some of these 
affects the outcome of our runs. Our models start with 
a common set of initial parameters summarized in Ta- 
ble n Most of these values are taken directly from ob- 
servations, but a few of the parameters deserve some 
additional discussion. Observations of the strength of 
magnetic fields in GMCs are quite uncertain. We set 
77b — 0.5, corresponding to equipartition between mag- 
netic and kinetic energy. For a more detailed analysis 
of the observati onal data and theoretical arguments for 
this choice, see iKrumholz fc McKed (^5, § 7.3). We 
set ?7v = 1.2 and (/)in = 1.0, corresponding to the de- 
cay of turbulence at the rate found by the simulations of 
iStone et"al] l)1998() . and to turbulence on the size scale of 
the entire GMC. Finally, we set 77E = 1, corresponding 
to HII regions injecting energy into turbulent media as 
efficiently as they would for smooth media. 

With these parameters fixed, we can fully specify 
the initial conditions for a model by giving the initial 
mass and column density of a cloud. We evolve mod- 
els with masses of Md.e — 0.2, 1.0 and 5.0, where 
Mci = Mc\fi X 10^ Mq- These masses span the range 
where most of the molecular g as in the Milky Way re- 
sides (|Wilhams fc McKeelll997li . We set the initial cloud 
column density to Nh,22 = 1-5, where Nh = Nh,22 x 10^^ 
cm^^, which is typical for Milky Wa y GMCs regardless 
of mass (.Larson 1981: Solomon et aUiigS? ) . We summa- 
rize the initial radius, velocity dispersion, and crossing 
time as a function of mass in Table El For reference, we 
also show how these quantities vary with column density 
at fixed mass. 

6. RESULTS 
6.1. Fiducial Runs 

We simulate the fiducial initial conditions 100 times 
each, using different random seeds for each run, for each 
of our three initial masses. Table |31 summarizes the sta- 
tistical outcome of our fiducial runs. 

The most basic result is that massive clouds attain a 
quasi-cquilibrium state, in which the decay of turbulence 
is roughly balanced by the injection of energy by HII 
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TABLE 2 
Initial cloud properties. 







Rci-o (pc) 


o-cl-o (km s 1) 


t„-o (Myr) 


0.2 


0.5 


33.7 


2.37 


13.9 


1.0 


0.5 


75.3 


3,54 


20.8 


5.0 


0.5 


168 


5.30 


31.1 


0.2 


1.5 


19.4 


3.12 


6.1 


1.0 


1.5 


43.5 


4.66 


9.1 


5.0 


1.5 


97.2 


6.97 


13.6 


0.2 


4.5 


11.2 


4.10 


2.7 


1.0 


4.5 


25.1 


6.14 


4.0 


5.0 


4.5 


56.1 


9.18 


6.0 



regions. In this state, cloud virial parameters fluctuate 
around unity, but most clouds spend most of their lives 
with virial parameters from 1 to 3, with a time-averaged 
value from 1.5 — 2.2 depending on the cloud mass, as 
shown by Figure ^ This may slightly overestimate the 
true virial parameter, because in our model we assume 
that all of the energy from HII regions goes into ran- 
dom turbulent motions that can then fuel cloud expan- 
sion, rather than into coherent motions of the cloud as a 
whole or, if HII regions fragment the cloud, into motions 
of these fragments. Such coherent motions are often ex- 
cluded in observational estimates of virial parameters. 
Nonetheless, at a qualitative level the result that clouds 
equilibrate to virial parameters ~ 1 is in good agree- 
ment with observations s howing that massive clouds are 
approximately virialized ijHever et al.l 120011 IBlitz et aP 

Hops). 

The evolution follows a sawtooth pattern, in which an 
HII region drives up the virial parameter, which then 
exponentially decays until it is increased by the next in- 
jection of energy. This is similar to the pattern seen in 
simulations by Li & Nakamura (2006), in which turbu- 
lence in star-forming clumps ~ 1 pc in size is maintained 
by energy injection from protostellar outflows. This equi- 
librium is maintained partly because the star formation 
rate responds to the current state of the cloud, increas- 
ing as the cloud contracts and its turbulence decays, and 
going back down when the cloud re-expands. Figure |21 
shows the depletion time, defined as the ratio of the cur- 
rent cloud mass to the current star formation rate, which 
exhibits this pattern. 

As shown in Figurc|31 cloud column densities also show 
a sawtooth pattern, oscillating up and down but remain- 
ing relatively constant over multiple expansion and con- 
traction cycles so that the time average is roughly the 
observed value Nh.22 ~ 1-5. (The figure is somewhat 
deceptive, because the longest lived-clouds tend to go to 
somewhat lower column densities, and these stand out 
the most when examining the figure because the shorter- 
lived clouds all overlie one another on the left side of the 
plot. The numerically-computed time-averaged density 
over all models is given in Tabled) Over their lifetimes, 
the average star formation efhciency of clouds, defined 
as the fraction of initial cloud mass transformed into 
stars, is 5 — 10%. HII regions ionize away anywhere from 
50 — 90% of the mass before finally destroying the clouds 
entirely, with lower mass clouds losing less mass to pho- 
toionization than more massive clouds. Figure ^ shows 
this evolution. 

The duration and stability of cloud equilibrium is af- 
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Fig. 1. — Virial parameter «vir versus physical time t and 
dimensionless time r for a sample of runs with fiducial parameters. 
We show Mci g = 0.2 (top panel), g = 1.0 (center panel), and 
Mj-i^g = 5.0 (bottom panel). 

fected by cloud mass. Low mass clouds, Mc\fi = 0.2 are 
only stable for an average of 1.6 crossing times (3.2 free- 
fall times), and are most often destroyed when they form 
an HII region that delivers enough momentum to com- 
pletely unbind them. They survive only one or a few cy- 
cles of expansion and contraction. In contrast, very mas- 
sive clouds, Mci_6 = 5, are stable for 3.2 crossing times 
and for multiple cycles before finally being unbound by 
HII regions. At all masses most clouds are destroyed 
by direct disruption rather than by a photodissociation, 
with the number photodissociated dropping sharply as 
a function of mass. This result is largely a function of 
cloud lifetimes. Direct disruption by a single HII region 
occurs when a truly large association forms, one that is 
on the tail of the mass distribution. Since larger clouds 
live longer and form more stars, they sample this tail 
more thoroughly, and thus are more likely to form one 
very large association capable of disrupting them than 
smaller clouds. An additional effect comes from the cloud 
velocity dispersion. Since larger clouds have higher ve- 
locity dispersions, the HII region shells driven by small 
associations merge more quickly, so the relative amount 
of energy injected by large versus small associations in- 
creases. 

Cloud disrup tion is much more fr e quent in o ur models 
than in those of lWilliams fc McKed 1)19971) or of lMatzneJ 
( 2002). The primary difference is a revised criterion for 
dynamical disruption. Williams & McKee took the on- 
set of a cometary phase to be the point at which an HII 
region effectively disrupts or displaces its parent cloud. 
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TABLE 3 
Fiducial run outcomes. 



Ma,6 


tiife (Myr) 


Ovir 


Nh,22 


SFE 




-^disrupt -^dissoc 




0.2 


1.6 (9.9) 


2.2 


1.4 


0.053 


0.59 


63 37 





1.0 


2.2 (20) 


2.1 


1.3 


0.054 


0.70 


92 8 





5.0 


3.2 (43) 


1.5 


1.5 


0.082 


0.80 


99 1 






Note. — Col. (2): Mean lifetime in crossing times (Myr). Cols. (3-4): Ovir and 
Nh,22 averaged over all times and runs. Col. (5): Star formation efSeiency. Col. 
(6): Fraction of mass photoevaporated prior to cloud destruction. Cols. (7-9): 
Number of runs out of 100 that ended in disruption, dissociation, and collapse. 




Fig. 2. — Depletion time versus physical time t and dimensionless 
time T for a sample of runs with fiducial parameters. We show 
M(.i g = 0.2 {top panel), M^i g = 1.0 {center panel), and Mj-i g = 
5.0 {bottom panel). 

Matzner considered the delivery of a momentum psh in 
excess of MdVcsc to define disruption. Our criterion 
(fsh > I'esc when r > Rd) resembles the Matzner thresh- 
old, but requires roughly half as much momentum be- 
cause we model HII regions as hemispher es and clouds as 
spheres. In addition, our adoptio n of the lKroupal 1)20021) 
IMF and lParravano et all l|2003D stellar properties leads 
to about twice as many ionizing photons per stellar mass 
relative to Williams & McKee and Matzner. The revised 
model of HII region expansion also plays a minor role. 
All told, the largest cloud that can be disrupted by HII 
regions is about ten times more massive in this work than 
in the Matzner models. Our disruption criterion, being 
more generous than the Williams & McKee or Matzner 
criteria, accounts in an approximate manner for cloud 
displacement and fragmentation by large HII regions. 
Our conclusion that more massive GMCs survive 



Fig. 3. — Column density Nh versus physical time t and di- 
mensionless time r for a sample of runs with fiducial parameters. 
We show M^i g = 0.2 {top panel), M^i g = 1.0 {center panel), and 
Mci^g = 5.0 {bottom panel). 

longer than les s massive ones also c ontrasts wit h the re- 
sults obtained iWilliams k McKed ^997) and iMatzneJ 
()2002(1 . One reason for our different result is that in our 
model the star formation rate per unit mass is propor- 
tional to the inverse of the free-fall time, which, at con- 
stant Nh, produces more star formation per unit mass in 
small clouds than in large ones. Thus, low mass clouds 
are subject to more vigorous feedback. In contrast, previ- 
ous studies generally adopted star formation laws under 
which the star formation rate per unit mass was inde- 
pendent of or increased with GMC size. This higher star 
formation rate compensates for the lack of large asso- 
ciations in small clouds and the ease with which small 
clouds can be driven into a cometary configuration, both 
of which tend to reduce mass loss. As a result clouds 
photoevaporate about the same fraction of their mass 
per dynamical time regardless of their starting mass. 
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Fig. 4. — Cloud mass {solid lines) and 10 times mass of 
stars formed {dashed lines) versus physical time t and dimension- 
less time T for a sample of runs with fiducial parameters. We 
show M^i g = 0.2 {top panel), A/j-i 5 = 1.0 {center panel), and 
M^l g = 5.0 {bottom panel). 

A second reason for our finding that small clouds live 
less time is our focus on the dynamical rather than 
the photevaporation lifetime. For massive clouds, our 
models show that the photoevaporation e-folding time 
is ~ 20 — 30 Myr, comparable to or a bit shorter than 
the dynamical lifetime, and consistent with earlier esti- 
mates. In contrast, for low mass clouds the dynamical 
destruction time is shorter than the photoevaporation 
time. Clouds are dynamically broken up by HII regions, 
either by direct disruption from a single region or ex- 
pansion to the point of photodissociation through the 
collective action of several HII regions, before most of 
their mass is photoevaporated away. This occurs because 
smaller clouds have less inertia and lower escape veloci- 
ties, so they are more easily pushed apart by HII regions. 
Larger clouds require truly large HII regions to unbind 
them, and so are not disrupted until they have lost sig- 
nificant mass and enough associations have formed to 
reach to this tail of the distribution. Smaller clouds can 
be destroyed by more typical HII regions. 

This distinction between photoevaporation and dy- 
namical lifetime points to an important caveat of our 
analysis. The dynamical lifetime we find is the time for 
which a cloud exists as a single coherent entity under 
the assumption of spherical symmetry. As we discuss 
above, it is possible that, rather than unbinding all the 
molecular gas, or expanding it to the point of photodisso- 
ciation, HII regions displace GMCs or blow GMCs into 
multiple pieces that are not bound to each other, but 



each of which remains molecular and internally bound, 
and can continue forming stars. Thus, the time for which 
the Lagrangian elements of a cloud remain molecular and 
star-forming may be longer than the dynamical lifetime 
of the cloud which we have found. The photoevaporation 
time may provide a better estimate of this Lagrangian 
lifetime. For clouds of mass Md.g ^ 1 this distinction 
is not that significant since the lifetime is comparable 
to the photoevaporation time, but it can be for smaller 
clouds. Determining whether parts of these clouds do 
continue forming stars even after they no longer form a 
single gravitationally bound entity will require radiation 
hydrodynamic simulations that can include both cloud 
dynamics and photoevaporation. 

6.2. Varying Column Density 

Prompted in part by the result that giant r nolec- 
ular clouds both in the Milky Way ijLarsonI Il981t 
ISolomon et al. 198 7f) and in nearb y galaxies with sim- 
ilar metallicitics ( Blitz et al.ll200!^ have a surprisingly 
small range of column densities, we investigate how vary- 
ing the initial column density affects the evolution. For 
each mass we run with the fiducial parameters, we re- 
run using the same 100 random seeds but initial col- 
umn densities of Nh.22 = 0.5 and 4.5. The low col- 
umn density case corresponds to a GMC with a mean 
column density that is just barely sufhcient to be self- 
shielding against the interstellar UV field. Giant clouds 
of such low column density (as opposed to isolated low- 
mass clouds) have not been seen in CO surveys either 
in the Milky Way or in other local group galaxies, and 
there is observational evidence and theoretical expecta- 
tion that molecular gas at column densities lo wer than 
Nn ^ 0.9 is not star-f orming (McKec 19 89t iLi et al 

■X . 1 . , 11 rTTTTTr? ^r\i^r\ nr^r^n t 1 i i i 



1997 : lOnishi et all '1998'. T999, 2002; J ohnstone et al 
20041 - though see Hatchcll ct al. 2005, who find less 



star formation at low column densities, but not a com- 
plete absence of star formation). However, in the model 
of GMCs as turbulent density fluctuations, low column 
density clouds of all masses are expected to exist and 
be star-forming (e.g. Vazquez-Scmadcni et al. 1997), so 
we consider the low column density case without impos- 
ing a column density threshold for star formation. The 
high column density case corresponds to GMCs found in 
galaxies that are entering the starburst regime. For ex- 
ample, in M64, a weak nuclear starburst galaxy, the typ- 
ical GMC column density is 2.5 times that in the Milky 
Way m^osolowskv fc R]it7il200l. 

Table 01 gives statistical results for the runs with vary- 
ing Nh, and Figure [S] shows the evolution of column 
density versus time for a sample of runs at each mass. 
The results show that Nh, 22 ~ 1 seems to be roughly a 
critical point in column density. Clouds that begin their 
evolution at substantially lower column density tend to 
disrupt or dissociate in ~ 1 dynamical time. Clouds that 
start at higher column densities show a pattern that de- 
pends on mass. At masses Afci,6 ^ 1, they remain stable 
for long times. At higher masses, they tend to undergo 
uncontrolled collapse. 

This change in behavior is not due to a change in 
the star formation rate per dynamical time. This is al- 
most independent of the column density, since SFRff oc 
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for fixed a^w- Instead, the evolution 
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TABLE 4 

Outcomes with varying column density. 







tlife (Myr) 




Nh,22 


SFE 


Afphot 


-^disrupt 


-^dissoc 


A'col 


0.2 


0.5 


0.44 (6.1) 


1.7 


0.60 


0.022 


0.34 


74 


26 





0.2 


1.5 


1.6 (9.9) 


2.2 


1.4 


0.053 


0.53 


63 


37 





0.2 


4.5 


5.6 (15) 


1.5 


5.2 


0.19 


0.80 


93 





7 


1.0 


0.5 


0.72 (15) 


1.9 


0.61 


0.026 


0.51 


8 


92 





1.0 


1.5 


2.2 (20) 


2.1 


1.3 


0.054 


0.70 


92 


8 





1.0 


4.5 
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83 


5.0 


0.5 


1.3 (41) 


1.6 


0.5 


0.039 


0.58 


5 


95 





5.0 


1.5 


3.2 (43) 


1.5 


1.5 


0.082 


0.80 


99 


1 





5.0 


4.5 


















100 



Note. — Column definitions are identical to those in Table l3l and cases with Nh,22 = 
1.5 are identical to the values in that Table. We compute average quantities excluding 
runs that result in collapse, and we do not attempt to compute averages in cases where a 
majority of runs produce collapse. 
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Fig. 5. — Column density A'^^f versus dimensionless time t for 
a sample of runs with varying initial Nh- Each panel shows a 
different initial mass, shown in the upper left corner. The dot- 
ted horizontal line indicates the column density at which clouds 
dissociate. 

significantly affects tlie efficiency of feedback;. Witli our 
approximation that merged HII regions do not inject en- 
ergy into clouds, the fraction of an HII region's energy 
that is lost to radiation increases with column density, so 
for equal ionizing fluxes and times, the kinetic energy in 

— 3/5 

an HII region shell varies as A'^^ ' . Furthermore, higher 
column densities increase the cloud velocity dispersion 
(as N^^^ for fixed avir) and decrease the expansion ve- 
locity, causing HII regions to break up and merge earlier. 
Since in our model HII regions only gain energy as long 
as they are expanding, the sooner they merge the less en- 
ergy they inject. From equation H35(l . for fixed cloud and 
HII region properties the ratio of the merger radius to 

the cloud radius scales as 7V^^^^°. Combining these two 
effects, equation (jHSfl shows that the energy injected by a 
single HII region of fixed properties into a cloud of fixed 
mass varies with column density as Nj^^^^^^^ . Thus, the 



efficiency with which ionizing photons are converted into 
turbulent motions is a reasonably strong function of the 
column density, and this picks out a particular character- 
istic column density at which energy injection balances 
loss. This turns out to be roughly the observed column 
density Nh,22 ~ 1-5. At this point we should mention 
one significant cautionary note: we have not explored 
the effects of the external environment of GMCs, and 
in particular the ambient pressure, on the characteristic 
column density. We will consider this effect in Paper II. 

It is easy to understand intuitively why the critical 
column density above which clouds tend to collapse for 

masses M^fi ^ 1 should be between Nh,22 — 1-5 and 4.5. 
In the high column density case, the velocity dispersion 
required to hold up a cloud is almost more than HII re- 
gions can provide. HII regions cannot effectively drive 
turbulence to velocity dispersions larger than the ion- 
ized gas sound speed cu = 9.7 km s^^, and as indicated 
in Table 13 in the high mass, high column density case 
the level of turbulence required to maintain avir = 1-1 
is (Jci-o = 9.2 km s^^. For the medium mass case it is 
(7ci-o — 7.0 km s^^. Furthermore, if the cloud contracts 
some so that its velocity dispersion increases a bit beyond 
these values, then HII regions will break up and merge 
with the turbulence almost immediatel y and will there- 
fore inject very little energy. Matzner ( 2002) previously 
discussed the inability of HII regions to maintain virial 
balance in systems where the required velocity dispersion 
exceeds cu, and our models find the same result. We dis- 
cuss the implications of this in more detail in 5 17.31 Note, 
however, that this conclusion is partially dependent on 
our approximation that HII regions cease driving turbu- 
lence once they cease dynamically expanding and merge, 
which may not be entirely correct - see § 17.41 

The long lifetime and high star formation efficiency 
produced at low masses and high column density also 
suggests an interesting interpretation: this combination 
of parameters may correspond to the regime of forma- 
tion of individual OB associations and open clusters, 
which is characterized by relatively high column den- 
sities l|McKee fc Tanll2003j) and star formation efRcien- 
cies > 10% llLa.da Kr. T.a,da l l200l . and probab ly requires 
many crossing times to complete l|Tan et al.l l2006'). Al- 
though the highest column density we have considered is 
still relatively modest by the standards of some cluster- 
forming clumps, which can reach Nh.22 of several tens. 
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the general result that low masses and high column den- 
sities can produce long-lived bound objects is suggestive. 
This is particularly true in light of the recent evidence 
that rich c lusters require ^ 5 crossing ti mes to assemble 
(iTan et al. 2006: Huff & Stahlcr 2006; K rumholz fcTal 
120061) ■ and must therefore be held up against collapse yet 
not be unbound by feedback for at least this long. 

6.3. Varying Dissipation Rate 

As we discuss in § 13.11 it is possible that the turbu- 
lent dissipation rate may vary from our fiducial esti- 
mate, either being substantially lower if simulations of 
turbulent decay have missed some important physics, or 
substantially higher if the characteristic scale of molecu- 
lar cloud turbulence is smaller than the size of an entire 
GMC. The latter possibility is probably ruled out by ob- 
servations showing that turbulence in molecular clouds 
is self-sim ilar out to the size of the entire GMC (e.g. 
lOsscnkop f fc Mac Lowll200l IHever fc BruntI 120041) . but 
we explore it nonetheless to better understand how the 
dissipation rate affects GMC evolution. We re-run our 
fiducial case rj^ — 1.2, (l)in = 1-0, corresponding to tur- 
bulence on the GMC scale and decay at the rate mea- 
sured by Stone et al. (1998), with rj^ = 0.4, (f>in = 1-0, 
corresponding to turbulence on the GMC scale decaying 
somewhat more slowly than Stone et al. find, and with 
r/v = 1.2, = 0.33, corresponding to turbulence driven 
on 1/3 of the GMC size scale decaying at the Stone et 
al. rate. Since the dissipation rate depends on the ratio 
Vv/f^in, we are therefore considering energy dissipation 
rates that are a factor of 3 smaller and larger than the 
fiducial case. (From eauations ll5l 1261 and 1281 this factor 
of 3 larger of smaller dissipation rate should correspond 
roughly to a factor of 3 change in the acceleration of the 
cloud radius, since £ (xij^a^, cr' oc £/cr^, and R" oc cr^.) 

Table [S] summarizes the statistical results of varying 
the dissipation rate, and Figure |S1 shows the evolution 
of column density versus time as a function of the dissi- 
pation rate and cloud mass. The general result is that, 
with the exception of the high mass, fast dissipation case, 
the results do not greatly depend on the dissipation rate. 
As the turbulent dissipation rate increases, the lifetime 
and mean virial parameter decrease and the mean col- 
umn density and the star formation efficiency increase, 
but none of these changes by more than ~ 10%. That 
changing the dissipation rate by a factor of three in either 
direction induces a much smaller change in the cloud col- 
umn density and virial parameter suggests that these val- 
ues represent roughly an equilibrium configuration, and 
that this equilibrium is quite robust. As the dissipation 
rate changes by a factor of nine from the lowest to the 
highest values we try, clouds contract a bit more, form 
stars a bit more vigorously, and are destroyed by stellar 
feedback a bit sooner, but only modest changes are suf- 
ficient to offset the change in dissipation rate. Star for- 
mation self-regulates to produce GMCs with iVH,22 ~ 1 
and ttvir ~ 1 — 2, and the regulation is stiff in the sense 
that even an order of magnitude change in the dissipa- 
tion rate does not alter it much. We discuss reasons for 
this in § O 

There are two exceptions to this, where varying the 
dissipation rate does produce a qualitative change in the 
outcome. For clouds of high mass and high dissipation 
rate, we see a phenomenon similar to the high mass, high 



column density case. For Afci,6 = 5, the velocity dis- 
persion starts at ctcI-o — 7.0 km s^^, but increases very 
rapidly as the cloud contracts due to decay of turbulence. 
By the time the first large HII regions have expanded to 
the point where they contain substantial kinetic energy, 
the cloud has contracted to the point where the veloc- 
ity dispersion required to hold it up is comparable to or 
larger than cn = 9.7 km s'-^. As a result, HII regions 
cannot hold up the cloud, and it collapses. 

The other exception occurs if we simultaneously in- 
crease the column density to Nh,22 = 4.5 and lower the 
dissipation rate by setting r^v = 0.4, = 1.0. In this 
case, clouds have long lifetimes and high star formation 
efficiencies regardless of their mass, and their column 
densities gradually decrease with time. Qualitatively, 
this is the same behavior we see for the case Mc\fi = 0.2, 
Nh,22 = 4.5, and the fiducial dissipation rate. Figure H 
shows the evolution of clouds within initial column densi- 
ties of Nh,22 = 4.5 as a function of mass and dissipation 
rate. As the plot shows, reducing the dissipation rate 
changes the characteristic mass at which one goes into 
the high star formation efficiency, long-lived regime. The 
boundary between this regime of very long stability and 
the regime of stability for a few dynamical times appears 
to depend weakly on both the mass and the dissipation 
rate. 

6.4. Varying Energy Injection Efficiency 

As we discuss in § 13.2.31 we are not certain of the effi- 
ciency with which HII regions are able to drive motions in 
turbulent media. We therefore re-run our fiducial models 
with r/E — 0.25, thereby reducing the amount of energy 
injected by HII regions by a factor of 4. This allows us to 
examine how strongly our results depend on our assumed 
efficiency. 

Table El summarizes our results, which show how small 
a difference the change in driving efficiency makes. With 
a factor of 4 less energy injection for an HII region of 
the same luminosity, the mean lifetime of the lower mass 
clouds we model increases by about 10% and that of 
the highest mass clouds decreases by the same fraction. 
Mean virial parameters decline by 10%, and mean col- 
umn densities rise a similar amount. Both star formation 
efficiencies and fractions of the mass photoevaporated 
rise, but again by only ~ 10%. The only quantity that 
changes significantly is the fraction of clouds destroyed 
by dissociation, which not surprisingly declines sharply. 
Thus, our results appear extremely insensitive to changes 
in the assumed energy injection efficiency of HII regions. 
Figure |H1 which shows the evolution of column density 
versus time for a sample of our runs, confirms this im- 
pression. The runs with lower energy injection efficiency 
go to slightly higher column densities, but overall show 
no major difference in their evolution from the fiducial 
case. Again, the star formation process seems to be self- 
regulating, so that large changes in efficiencies, either 
of energy injection or of dissipation, produce only very 
small changes in the global statistics of cloud evolution. 
We discuss reasons for this in § 17.11 

7. DISCUSSION 

7.1. The Relationship of Cloud Properties to Rates of 
Radiative Loss and Energy Injection 
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TABLE 5 

Outcomes with varying dissipation rate. 



,6 (pin tlifo (Myr) cEvir Nh,22 SFE Afphot A^dis rupt -^dissoc -^col 



0.2 


0.4 


1.0 


1.7 (10) 


2.8 


1.2 


0.046 


0.53 


49 


51 





0.2 


1.2 


1.0 


1.6 (9.9) 


2.2 


1.4 


0.053 


0.59 


63 


37 





0.2 


1.2 


0.33 


1.4 (8.4) 


1.9 


1.8 


0.067 


0.57 


44 


56 





1.0 


0.4 


1.0 


2.3 (21) 


2.6 


1.1 


0.046 


0.64 


90 


10 





1.0 


1.2 


1.0 


2.2 (20) 


2.1 


1.3 


0.054 


0.70 


92 


8 





1.0 


1.2 


0.33 


1.7 (16) 


1.7 


2.1 


0.074 


0.74 


61 


35 


4 


5.0 


0.4 


1.0 


4.4 (60) 


2.3 


1.0 


0.070 


0.75 


98 


2 





5.0 


1.2 


1.0 


3.2 (43) 


1.5 


1.5 


0.082 


0.80 


99 


1 





5.0 


1.2 


0.33 


















100 



Note. — Column definitions are identical to those in Table O and cases with r)v = 1-2, 
(/)in = 1.0 are identical to the values in that Table. As in Tablel4l we compute average quantities 
excluding runs that result in collapse, and we do not attempt to compute averages in cases 
where a majority of runs produce collapse. 



E 



E 




Fig. 6. — Column density versus dimensionless time t for a sample of runs with varying dissipation rates and masses. The mass and 
dissipation rate are indicated in each panel - slow is rjv = 0.4, 0in = 1.0, medium is rjv = 1-2, = 1.0, and fast is rjv = 1.2, tpi^ = 0.33. 
The dotted horizontal line indicates the column density at which clouds dissociate. All runs begin with Nh,22 = 1-5. 
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Fig. 7. — Column density Nu versus dimensionless time r for a sample of runs with varying dissipation rates and masses, all starting 
from high initial column densities. The mass and dissipation rate are indicated in each panel - slow is rjv = 0.4, = 1.0, medium is 
7?v = 1.2, = 1.0, and fast is r)v = 1.2, (f>\^ = 0.33. The dotted horizontal line indicates the column density at which clouds dissociate. 
AH runs begin with Njj 22 = 4.5. 



TABLE 6 

Outcomes with varying HII driving efficiency. 







tlifc (Myr) 


Qvir 


Nh,22 


SFE 


A^phot 


-^disrupt 


-^dissoc 


A'col 


0.2 


0.25 


1.7 (11) 


1.6 


1.7 


0.065 


0.67 


93 


7 





0.2 


1.0 


1.6 (9.9) 


2.2 


1.4 


0.053 


0.59 


63 


37 





1.0 


0.25 


2.1 (20) 


1.6 


1.6 


0.062 


0.75 


100 








1.0 


1.0 


2.2 (20) 


2.1 


1.3 


0.054 


0.70 


92 


8 





5.0 


0.25 


2.8 (39) 


1.3 


1.8 


0.087 


0.82 


100 








5.0 


1.0 


3.2 (43) 


1.5 


1.5 


0.082 


0.80 


99 


1 






Note. — Column definitions are identical to those in Tablel^ and cases with r/E = 1.0 
are identical to the values in that Table. 
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Fig. 8. — Column density Nh versus dimensionless time r for a sample of runs with varying HII region driving efficiencies and masses. 
The mass and driving efficiency are indicated in each panel. The dotted horizontal line indicates the column density at which clouds 
dissociate. All runs begin with Nu 22 = 1-5. 



One of the most interesting results of our analysis is 
how little the globabl statistics of cloud evolution, such 
as lifetimes and star formation efhciencies, change in re- 
sponse to large changes in parameters such as the rate 
at which turbulence decays via isothermal shocks or the 
efficiency with which HII regions transfer their kinetic 
energy into turbulent motions. We can understand intu- 
itively why large variations in either driving efficiency or 
dissipation rate cause such minimal changes in the stati- 
tistics of cloud evolution by examining Figure |21 As the 
Figure illustrates, clouds form most of their stars during 
their contraction phases, when they reach high densities, 
so the rate of star formation is effectively set by the ca- 
dence of expansion and contraction cycles. This cadence 
is affected only slightly by the dissipation rate and the 
energy injection efficiency, because both of these prop- 
erties are coupled to cloud evolution poorly, with large 
time delays. 

For a GMC to contract from an expanded state, the 
rate-limiting step is not the time required to dissipate 
the energy injected by HII regions, it is the cloud free-fall 
time, the time the cloud would take to collapse even if it 
contained no turbulence at maximum expansion. Thus, 
varying the turbulent dissipation rate makes little differ- 
ence to the cloud contraction time. Similarly, the time it 
takes for HII regions to drive apart a cloud is controlled 
less by the amount of energy added per unit mass of stars 
formed than by the time delay between when stars begin 
forming and when the HII regions they create expand, 
break up, and drive turbulent motions. Lowering the 
driving efhciency means that clouds expand somewhat 
less far, but since the dissipation rate of the turblence is 



a strong function of the velocity dispersion, £ cx cr'^, the 
extra energy at early stages produced by higher efficiency 
is radiated away very quickly, and produces only slightly 
increased cloud expansion. We can put this argument 
more formally by noting that our energy equation gives 
a' oc (7^ , neglecting changes in velocity dispersion due to 
external pressure and cloud expansion, so the time re- 
quired for a velocity dispersion of ctq immediately after a 
large HII region breaks up to decay to a value ai obeys 
t cx (fJo — fij/lcoci). If the final velocity dispersion is 
much smaller than the initial one, i.e. ui ctq, this ratio 
is independent of ctq. Thus, the decay time for a large 
amount of initial energy injected is only slightly larger 
than the decay time for a significantly smaller energy in- 
jection, so the cloud expansion time and maximum radius 
depend only weakly on the driving efficiency. 

This insensitivity to the dissipation rate and the en- 
ergy injection rate breaks down if the dissipation rate 
becomes so high or the energy injection so inefficient that 
HII regions are not capable of expanding clouds at all. 
In this case, clouds undergo runaway collapse. It also 
breaks down if the dissipation rate becomes so low or 
the efficiency so high that the first generation of HII re- 
gions simply unbinds most clouds, in which case clouds 
form only one generation of stars and are destroyed on 
a time scale set by the expansion time for HII regions 
from that first generation. However, in between these 
two extremes there is a broad range of parameter space 
within which the cadence of expansion and contraction 
cycles is set mostly by the cloud free-fall time and the 
time required for HII regions to break up, with only a 
very weak dependence on the details of the energy loss 
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rate or the energy gain efficiency. The degree of insen- 
sitivity is illustrated by Figures El and |H| in which factor 
of several variations in dissipation rate or driving effi- 
ciency change the expansion and contraction period by 
only tens of percent. 

7.2. GMC Stability and Lifetime 

Observations place strong constraints on the stability 
of giant molecular clouds, and are now beginning to con- 
strain their lifetimes as well. We must test theories that 
seek to explain the behavior of GMCs against these ob- 
servations. The first constraint is that giant molecular 
clouds cannot be in a state of global collapse. If they 
were, then they would convert order unity of their mass 
into stars in a crossing time, and this would produce a 
star formation rat e two orders of magnitude l arger than 
the observed one ijZuckerman fc Evan^ll974j) . Thus, a 
model of GMC evolution must explain why GMCs are 
stable against global collapse and convert only a few per- 
cent of their mass into stars per crossing time. 

The second observational constraint is that GMCs, at 
least the most massive ones, live for considerably more 
than a single crossing time. While GMC lifetimes are 
quite difficult to estimate inside the Milky Way, age 
spreads of the largest OB associations are ~ 20 Myr 
l)Blaauwl lT964: Bhtz & Shu 1980), which is a probable 
lower limit on the lifetimes of GMCs. More robust con- 
straints are available in extragalactic observations, where 
there is no distance ambiguity. Associations between 
GMCs and star clusters imply a typical GMC lifetime 
of ~ 20 Myr i n M33 (Enaarriola et al. 2003) and 27 Myr 
in the LMC f Fukui et a l. 1999; Blitz ct al. 2005). This 
is significantly greater than the GMC crossing time of 

^ 10 Myr. While the difference between one and a few 
crossing times might seem insignificant, recall that the 
e-folding time for the decay of turbulence is roughly a 
crossing time. A cloud two crossing times old will have 
lost almost all of its turbulence if feedback or some other 
source cannot replenish it, and will therefore undergo 
global collapse. Thus, the observations not only require 
that clouds be stable against collapse, they require that 
this stability be maintained for several turbulent decay 
times. 

The GMC model we present in this Paper, using the 
fiducial parameters suggested by observation and pre- 
vious theoretical work, provides very good qualitative 
agreement with the observations. We show that feed- 
back from star formation can keep clouds supersonically 
turbulent and virialized for ~ 30 Myr, until they are 
destroyed by feedback. Combined with the results of 
iKrumholz & McKec ( 2005) showing that supersonic tur- 
bulent motions naturally produce a star formation rate of 
a few percent per crossing time if star formation occurs in 
virialized clouds, this model satisfies both observational 
constraints. We discuss the question how how GMCs 
evolve during this lifetime in more detail in § 17.31 

Alternative models run into difficulty with one of these 
two observations, or with other data. One possible ex- 
planation why GMCs do not undergo global collapse is 
that they are gravitationally unbound transient fluctu- 
ations in the atomic ISM; only local subregions consti- 
tuting a sma ll fraction of the clo ud mass are bound and 
can coUpase llMac Low fc KlessenL2004:. Clark fc Bonnelll 
l2?ffli IClark et all 120051 IVazauez-Semadeni et all 120061: 



IDobbs et al.ll200l . However, in this case it is hard to 
see how GMCs cou ld survive 20 — 30 Myr. For example, 
IClark et alJ |2005l) find GMC lifetimes of only about 10 
Myr in their simulations of unbound clouds. There are 
also other severe observational difficulties. There are no 
molecular clouds with masses ^ 10* Mq that are clearly 
gravitationally unbound, wi th virial parameters c tyir 3> 
1, either in the Milky Wav iHeyer fc Brun t"2004» or in 
other galaxies fe ngargio l^ et alJl2003nRo solowsk v et alJ 
'2003 URosolowskv fc Blit z 2005; Blitz ct al. 2005)." There 
are many examples in the Milky Way of clouds of mass < 
IO^Mq with virial parameters avir 3> 1 ( Hevcr fc Brui3 
t2-QQ4t) , and there is no obvious reason why massive clouds 
should not also display a range of avir if they are gen- 
erally unbound. Transient fluctuation models also have 
problems explaining the existence of a GMC mass scale 
and the low rotation rat es of GMCs. These argume nts 
are presented in detail in IKrumholz fc McKed l|2005fl . 

A second possiblity is that GMCs are bound, at least 
marginally, but that they are destroyed before the tur- 
bulence with which they are born decays away. As a 
result, they never have a chance to begin global collapse. 
However, the observed lifetimes of ?J 20 Myr appear to 
rule out the possibility that such rapid destruction is 
the norm. Furthermore, to work this model requires a 
mechanism of cloud destruction that reliably operates 

in ^ 1 crossing time. As we show here, HII regions 
arc not capable of completely photoevaporating or dis- 
rupting massive clouds over such short periods. Super- 
novae and protostellar winds inject consierably less en- 
ergy into G MCs per unit mass of stars formed than do 
HII regions llTenorio-Tagle et al.llT985t lYorke et alllTgSgt 
lMatzn er"200^. so they are unlikely candidates for rapid 
cloud disruption. In the absence of a plausible mecha- 
nism for cloud disruption in ^ 10 Myr, this model faces 
a major theoretical problem as well as an observational 
one. 

It is worth noting at this point that arguments for 
GMC lifetimes in the Milky Way strict ly limited to 1 
crossing time, e.g. iHartmann et alJ l)2001j) . generally rely 
on observations of clouds in the solar neighborhood, all 
of which are much smaller than the larger, more typical 
GMCs we have considered in this paper. Since we find 
that molecular clouds ^ 10^ Mq in mass do only sur- 
vive for ^ 1 crossing time, our results are consistent with 
the idea that solar neighborhood GMCs are short-lived. 
We simply suggest that GMC lifetime is mass-dependent, 
and that the nearest clouds, due to their atypical masses, 
also have atypical lifetimes. 

7.3. An Evolutionary Scenario for GMCs 

Our findings allow us to present a rough evolutionary 
scenario for GMCs that is consistent both with the con- 
straints of stability and GMC age described in S 17. 21 and 
with the observation that GMCs in the Milky Way and in 
other galaxies all lie in a narrow range of column densi- 
ties and virial param eters, ce nterin g around Nh.?.? ~ 1.5 
and avir » 1 - 2. l|Larsonlll98lt ISolomon et all 119871: 
IBlitz et al.ll200l. T his scenario is quite similar to that 
suggested in lMcKe^ (|1989ri and is developed in consider- 
ably greater detail in Paper II. 

Our scenario is that GMCs are born at low column 
densities by condensation out of the atomic ISM, pri- 
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marily triggered by self-Kravitational instabilities in spi- 
ral sh ocks fe.g. iKim fc Ostrikeil l20?nt iKim et all 12001 
l2003jl . This can only occur in regions of a galactic disk 
that are at significantly higher densities and pressures 
than is typical of the ISM. If star formation were able to 
start in such clouds immediately, they would undergo 
rapid disruption, converting ^ 3% of their mass into 
stars. However, in reality such clouds are probably not 
star-forming over most of their volume, so feedback is 
suppressed. 

This may explain the observation that GMCs pass 
through a phase in which they appear to lack embed- 
ded HII regions. The duration of this phase is diffi- 
cult to obtain fr o m ob servati ons - J ukui et al. (1999'), 
lEngargiola et alJ l|200.'^ and iBlitz et all l|200a[ ) report 
that ~ 1/4 of GMCs do not show Ha signatures asso- 
ciated with HII regions, but this is probably an overes- 
timate because it does not include highly obscured HII 
regions. Although the Ha sensitivities of the catalogs 
used by these authors are sufficient to detect HII regions 
with Ha luminosities comparable to the Orion Nebula, 
Orion is visible largely because it is on the near side 
of a dark cloud; at the distance of the LMC or M33, 
Orion would be invisible if it were oriented with the dark 
cloud on the near side rather than the other way around. 
Spitzer/MIPS observations have failed to find any GMCs 
that are dark at 24 /im, so there must be some embed- 
ded star formation even in GMCs without detected HII 
regions (E. Rosolowsky, 2006, private communication). 
This result suggests either that the HII regions surveys 
are incomplete, or that clouds do not begin making mas- 
sive stars until well after they have begun forming low 
mass stars, an effect we have not considered. 

Regardless of how long this initial starless phase lasts, 
a question we will address using our models in Paper 
II, in the absence of feedback, GMCs will contract due 
to loss of turbulent support, gravity, and external pres- 
sure, until they approach Nh,22 ~ 1- At that point, 
they become star-forming, and feedback stabilizes them 
against further contraction and keeps them in virial equi- 
librium for several crossing times. We observe GMCs at 
Nh,22 ~ 1, ttvir ~ 1 — 2 because this is where they spend 
the vast majority of their lifetimes. This column density 
is selected because it is the one for which energy injection 
by turbulence roughly balances energy loss by isothermal 
shocks. In addition to turbulent energy injection, the re- 
coil momentum produced by mass being evaporated by 
HII regions may also play a significant role in confining 
clouds. 

This quasi-equilibrium endures for an amount of time 
that depends on the cloud mass. For massive clouds, 
which contain most of the molecular mass in the Milky 
Way and in all but one other galaxy in which it is possible 
to estimate cloud mass distributions iFukui et al.llT999l: 
lEngargiol a eral]l2003t iBlitz et alJl2005jr it endures 2 — 3 
crossing times, or 20 — 30 Myr. During a cloud's life- 
time, it converts 5 — 10% of its mass into stars. These 
results are robust against changes in the assumed dissi- 
pation rate for turbulence or efhciency of turbulent driv- 
ing by HII regions. Less massive clouds are probably 
dynamically disrupted in ~ 1 crossing time, consistent 
with estimates of the l ifetimes of small clouds in the so- 
lar neighborhood fe.g. iHartmann et al.ll200'll) . although 



parts of them may endure in molecular phase and con- 
tinue forming stars for longer periods of time even after 
the original cloud has been disrupted. 

One interesting question is how our results might 
change if we considered a galaxy quite different from the 
Milky Way. Observations indicate that GMCs in normal 
spiral galaxies like the Milky Wa y are generally qui te sim- 
ilar to those in the Milky Way ijBlitz et al.ll2005|) . so to 
find truly different conditions we must consider galaxies 
that are approaching the regime of starbursts. The sole 
observational example of such a galaxy where we have ob- 
served the clouds is M64, a weak starburst in which the 

largest GMCs are ^ lO'' Mq in mass, and have surface 
densities up to Nh ~ 4 x 10 cm^^ (Rosolowsky & Blit3 
12005(1 . Despite these differences from Milky Way GMCs, 
observations indicate that these clouds are still in ap- 
proximate virial balance, and that like Milky Way GMCs 
their star formation rate is only a few percent per free-fall 
time. Thus, they are not in a state of global collapse. 

Since we have found that HII regions cannot sus- 
tain the turbulence and prevent collapse in such clouds 
(though see S I7.4|I . we are left with two possibihties. Ei- 
ther the dissipation rate in such clouds is lower than our 
fiducial estimate, or there is an additional source of en- 
ergy input that supplements HII regions. One strong 
candidate for a source of ene rgy injection is drivin g of 
turbulence by external shocks l|Kornreich fc Scalol2000 '). 
either from supernovae, gravitational instabilities driven 
by the po tential of the stars, or cloud-cloud collisions 
(Tan 2000). This mechanism encounters considerably 
difficulty in the Milky Way because GMCs are much 
denser than the gas around them in which shocks pro- 
pogate. This creates an "impedance mismatch" that 
makes it difficult to d rive turbulence into the GMCs 
(|Nakamura et al. 2006). However, in a starburst where 
the ISM is entirely molec ular, clouds are not much d enser 
than their surroundings. iRosolowskv fc Blitd l|2005fl find 
that in M64 the GMCs are overdense only by factors 
^ 2. This removes the impedance mismatch problem, 
and makes it far easier for external shocks to drive tur- 
bulence than it is in galaxies like Milky Way. Whether 
this mechansim can work in detail will require more an- 
alytic models and simulations to determine. 

7.4. Limitations of This Model 

The most obvious limitation of our model is the con- 
straint that GMCs evolve homologously. In mathemati- 
cal terms, this constraint is equivalent to dropping time 
derivatives of aj in the virial theorem. This should only 
change our results substantially if changes in the moment 
of inertia of a CMC occured primarily through changes 
in its shape rather than overall expansion or contrac- 
tion. While a cloud in approximate equilibrium probably 
does experience most changes in its inertia via changes 
in shape, a cloud that is undergoing global collapse or 
global disruption would almost certainly experience most 
changes in its inertia through those processes rather than 
through changes in shape. As a result, our broad conclu- 
sion that neither overall collapse or disruption occur for 
several crossing times seems likely to be robust. 

A more subtle limitation to our model is our simple 
boundary conditions for GMCs. For simplicity we have 
assumed a fixed mass budget and a fixed external pres- 
sure that puts GMCs into pressure balance initially. In 
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reality, GMCs may start forming stars while still ac- 
cumulating gas. The somewhat higher mean mass for 
GMCs with associated HII regions tha n for GM Cs with- 
out such associ ations seen in the LMC ijFukui e t al. 19995 
iBhtz et alf m)5) seems to point in this direction. In ad- 
dition, GMCs form in spiral shocks. As a result, the 
Lagrangian mass elements from which a cloud forms are 
probably subjected to a rising pressure as they enter the 
arm region, and then experience a falling pressure as the 
cloud moves out of the arm. This may have important 
effects on GMC evolution, and we explore the effects of 
more complex boundary conditions in Paper II. A related 
issue in GMC evolution, which we will also explore, is the 
effect of the apparent existence of column density thresh- 
olds for star formation. If clouds cannot begin forming 
stars until they accumulate a certain column of molecu- 
lar mass, then feedback in GMCs will not start up until 
some time after the cloud first becomes molecular. This 
lag may affect GMC evolution. 

There are also several limitations associated with our 
treatment of HII regions. We assume that the mass of the 
cloud is large enough that its evolution is dominated by 
HII regions, not protostellar outflows. We assume that 
the clouds are spherical, which means that all the en- 
ergy injected by HII regions goes into internal motions, 
whereas in fact some of it goes into moving the entire 
cloud. Finally, we neglect possible energy injection by 
HII regions after merging. In our model, shells driven 
by small HII regions or those in very turbulent clouds 
may drop to expansion velocities below the cloud veloc- 
ity dispersion well before the driving stars burn out. At 
this point the HII region will cease to drive a coherent 
expanding shell, and we approximate that the energy in- 
jection ceases. However, if the driving stars continue 
to ionize a part of the cloud, there is another potential 
driving mechanism. Clumps of gas that enter the ionized 
region will be photoevaporated on one side, and there- 
fore will be rocketed away from the ionizing stars. Even 
though the resulting thrust is not enough to drive the 
clumps out of the HII region entirely, since the expan- 
sion velocity is smaller than the cloud velocity dispersion, 
it can still alter the tr ajectories of gas clumps and po- 
tentially inject energy l|Tanll2r)ni^ . This effect could be 
important in the cases where we find that HII regions 
cannot support GMCs because the velocity dispersion 
required to hold up the cloud is comparable to or larger 
than the ionized gas sound speed. Determining whether 
this effect is significant will require either more detailed 
theoretical treatment or numerical simulations. Despite 
these caveat s, th ough, our results in § 16.316.41 and our 
analysis in S IT.ll show that our conclusions regarding the 
lifetime and evolutionary history of GMCs are quite ro- 
bust against changes in the assumed efficiency with which 
HII regions drive turbulent motions or the rate at which 
those turbulent motions decay. 

8. CONCLUSIONS 

In this Paper we derive the basic equations governing 
the evolution of the mass, radius, and velocity disper- 
sion of giant molecular clouds under the approximation 
of homologous motion. We construct simple models for 
the rate at which turbulent motions decay due to radi- 
ation from isothermal shocks and for the rate at which 
HII regions driven by massive stars within the cloud drive 



turbulent motions, and we use these to study the global 
evolution and energetic balance of GMCs. This enables 
us to build GMC models in which we neither assume 
energetic and virial balance, nor neglect the effects of 
feedback. Thus, we are able for the first time to address 
critical questions such as whether GMCs are in virial or 
energetic equilibrium, how long they live, what deter- 
mines their lifetimes and star formation histories, and 
what determines properties such as column density and 
virial parameter that appear to be roughly the same for 
aU GMCs. 

Our primary conclusion is that giant molecular clouds 
with observed properties are indeed quasi-equilibrium 
objects. The rate at which they lose energy via radiation 
from isothermal shocks is roughly balanced by the rate 
at which HII regions from star formation inject energy 
back into the cloud. This feedback keeps them virialized, 
ctvir = 1 — 2, and keeps them at the observed column den- 
sity Nh ~ 1.5 X 10^^ cm~^. Our results suggest that this 
column density, which appears to be generic for GMCs 
both in the Milky Way and in other galaxies where GMCs 
are observable, is in fact the result of feedback, since the 
efficiency of star formation feedback varies with column 
density, and the observed column density corresponds 
to the one required for equilibrium. Whether the GMC 
formation process, particularly the passage through an 
overdense and overpressured spiral shock, also plays a 
role in selecting the GMC column density is at this point 
an open question, one we will address in Paper II. 

The duration of the equilibrium state of a GMC de- 
pends on cloud mass, but for clouds > 10^ M© m mass 
it lasts for 2 — 3 crossing times, or roughly 30 Myr. The 
dominant destruction mechanism for GMCs is dynamical 
unbinding by the momentum delivered by an expanding 
HII region, but for massive clouds this unbinding does 
not happen until HII regions have photoionized away 
~ 90% of the GMC mass. Smaller clouds are dynam- 
ically disrupted by HII regions in ^ 10 Myr while ~ 50% 
of their mass is still in molecular form, but it is unclear 
if this disruption actually destroys their molecular mate- 
rial or simply breaks it into smaller, mutually unbound 
clouds. 

Overall, our models produce a picture of GMCs that 
is quite different than that suggested by models in which 
GMCs are gravitationally unbound or in which feedback 
from star formation is neglected. Models that neglect 
feedback appear incapable of reproducing the observed 
lifetimes and properties of the GMC population. We 
conclude that any reasonable model of GMC evolution 
cannot neglect the effects of star formation feedback, and 
that even a simple model of feedback such as the one 
we use here produces good agreement with the observed 
properties of GMCs. 
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APPENDIX 

A. DERIVATION OF THE VIRIAL THEOREM FOR AN EVAPORATING CLOUD 

Here we derive the e quation of motion for an evaporating homologous cloud, which is a generalization of the Eulcrian 
Virial Theorem (EVT. 'McK ee fc Zweibell992(l . Gas in the cloud is injected into the wind at a rate p, so the continuity 
equation for the cloud is 

|^ = -V.pv + p, (Al) 

where p and v are understood to be the density and velocity of cloud and not wind material. The wind satisfies its 
own continuity equation dpw/dt = —V • Pw^w — P, so that the summed equation d{p + pw) — —V • (pv + Pw^w) is 
simply the ordinary continuity equation. The equation of momentum conservation for the cloud gas is 

^(pv) = -V • (n - Tm) + pg + p(v + v^j), (A2) 

where 11 is the gas pressure tensor, Tm = [BB — (l/2)i?^I]/(47r) is the Maxwell stress tensor, B is the magnetic field, 
and g is the gravitational force per unit mass. Intuitively, the pv term represents the change in cloud momentum due 
to mass being transferred into the wind, and pv^ is the recoil momentum from the ejection. As with the continuity 
equation, there i s a corresponding m omentum equation for the wind. 

We now follow iMcKee fc Zweibell l(l992,') in deriving the EVT for a cloud with a wind. The moment of inertia of the 
cloud is 

/ci= / pr^dV, (A3) 

where the fixed volume of integration V^ir is chosen to be sufficiently large that it includes the entire cloud at all times. 
Taking the time derivative 

id = - / (V • pv)r2 dV + j pr"^ dV 

= - I {pvr^) ■dS + 2 I pvrdV + aiAURl^ (A4) 



where S'vir is the surface of the volume V^w, and 

nj = 

5-fc 



a, . ^A. (A5) 



Differentiating again, 

+ ^aiMcii?ci + fliMcii?ci4i. (A6) 

The time derivative can be taken out of the integral in the first term because S'vir is fixed, and for the second term 
we can substitute for the integrand using the momentum equation (|A2I) . This gives 

^/ci--^^ X '''^^ ' " X • ■ " - Pv • g} dV 

+ [ pr-vdV+ [ prv'^jdV + l-aiAURli + aiMciRciRci- (A7) 

To evaluate the term Jy pr • v dV, we make use of our homology assumption, which allows us to write the velocity 
at any point in the cloud as 

v^r^^+vturb, (A8) 

-flcl 

where the first term represents homologous expansion or contraction of the cloud, and the second is a turbulent velocity 
that carries no net radial flux of matter. Thus, /y . pr • Vturb dV = 0, and /y . pr • v = Jy pr^Rd/Rd- 
This allows us to write the final EVT, 



Svi, 

2aiA/cii?ci4i + |— |^Mcii?ci«;j + ^aiMcii?ci- (A9) 
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In this equation, the kinetic term is 

r = 

2 



\ I {3P + pv')dV (AlO) 



Kir 



To = i / r U dS, (All) 



for gas thermal pressure P, the surface term is 
the gravitational term is 

[ pr-gdV, (A12) 

and the magnetic term is 

^-^/^ (s2-Bo')rfV^, (A13) 

where Bq is the background magnetic field far from the cloud, and we require that Vvir be large enough to include the 
full volume over which the background field is distorted by the presence of the cloud. 

B. DERIVATION OF ENERGY CONSERVATION FOR AN EVAPORATING CLOUD 

Here we derive the equation of energy conservation for an evaporating homologous cloud. First, we rewrite the 
momentum equation (jA2|) in a slightly different form, 

d-v J X B , , , 

p— = -WP^pWcl>+^—+pv'^^, (Bl) 

where the terms on the right hand side are the pressure, gravitational, and Lorentz forces, (j) is the gravitational 
potential, J is the current density, and we have replaced the pressure tens or 11 with the isotropic pressure PI under 
the assumption that viscosity is negligible. Taking the dot product of (|B1|I with v yields 



§1 {If"'') + ^ ■ {l''''''^) = -V ■ VP - pv V0 + ^ • (J X B) + p + V • . (B2) 

Using the continuity equation, we can rewrite the gravitational work term as 

- pv • V0 = - V • pv(t> + <j)\7 ■ pv (B3) 

= - V • pv(/) - - +P^+ (B4) 

Similarly, we can rewrite the magnetic work term using Poynting's Theorem, which in the MHD approximation is 

^^(P2 - ^ ^ . ^ _j . E = -X . (J X B), (B5) 

where Sp is the Poynting flux, and we have used the fact that Bo is constant to include it in the time derivative for 
future convenience. Substituting into (|B2|I gives the equation for the time-evolution of the non-thermal energy, 

|Qp.2+p0+^?^^)+V.pvQ.2 + <^)+V.S, = -v.VP + p^+pQ.2 + + v.v^). (B6) 

To include the internal energy, we write down the first law of thermodynamics, 

de 

p- + PV-v = T-A, (B7) 

where e is the internal energy per unit mass and F and A are the rates of radiative energy gain and loss per unit 
volume. Combining this with the continuity equation gives the evolution equation for the internal energy, 

d / P\ 

— (pe) + V • pv f e + - j = pe + V VP + F - A. (B8) 

Adding together the evolution equations (jB6p and ljB8|) for the non-thermal and thermal energies gives the total 
energy equation 



d_ 

di 



2 Stt 

1 



V • pv ( + e + ^ + ) + V • Sp 



P^+p[^v^ + e + (t> + ^-^',;]+T-k. (B9) 
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To derive the global form of the energy equation, we integrate over the virial volume Vvir, which gives 



Itt ~ 

_ 1 
~ 2 



1 



-{£-B)- 



dS 



Sp • dS 



A/cli?clWoj 



where 



£ 



Stt 



(BIO) 



(Bll) 



is the total energy in the cloud, Qc\ and Cc\ are the total rates of radiative energy gain and loss integrated over the 
cloud, and we have used our assumption of homologous motion to evaluate the terms involving p and Vgj . To evaluate 
the first integral on the left-hand side, we note all the terms proportional to density vanish at the cloud surface because 
the ambient material has zero density, but that constant pressure and zero density correspond to an incompressible 
fluid, so the velocity of the ambient medium across the virial surface is related to the velocity of the cloud surface by 



i?ci(-Rci/i?vir)- This implies that 



Svi, 



PI 2- 



P 



V dS = 47rFambi?cl-Rcl. 



(B12) 



To evaluate the integral on the right-hand side, recall that gas in the wind contributes negligibly to the gravitational 
potential. Since the potential arises solely from cloud material, we can write (.,Shu. l99^ 



Vvir 



dt 



2 dt 



>dV 



Kir 



Mr- 



-W. 



(B13) 



Finally, we must evaluate Sp ■ dS, which represents the flux of magnetic energy across the virial surface as 
it is carried off by the wind. This is uncertain, because it depends on the process by which the mass is removed. 
As in S 12.11 we write the magnetic energy as a sum of turbulent and non-turbulent contributions, l/(87r) Jy (B^ — 

Bq) dV = B = Bnon-turb + 'Bturb- For the turbulcnt part, we assume that the wind does not change the nature 
of the MHD turbulence within the cloud, so the turbulent magnetic energy remains proportional to the turbulent 
kinetic energy. Since the loss of turbulent energy accompanying mass loss is Jy 2>pa'^/2dV = (Mci/A^ci)'?turb, this 

requires that the flux of turbulent magnetic energy be dBtmh/dt = (Mci/Mci)Sturb. For the noi i-turbulent part , for 
simplicity and based on the observation that GMCs have roughly constant mass-to-flux ratios l|Crutcherlll999D . we 
assume that the mass-to- flux ratio remains constant as the cloud loses mass. Since Bnon-turb oc this implies 
dBnou-tmh/dt = 2($/<I>)Bnon-turb = 2(Afci/Mci)Snon-turb- Thus, wc Can Write 



Sp • dS = -^{Btn^b + 2Bnon-turb) = ~^{B ~ T^IW). 
Mc\ Mc\ 



Substituting l|B12|l . I|B13(I . and l|B14|l into (|B10(I . we arrive at the flnal energy equation for the cloud: 

k \ ■ ■ 

Mcii?ci<i+eci-/:ci. 



'^^ [£+{1- vim - 4^^'a.„bi??i4l + 



dt Mc 



4 -A:, 



(B14) 



(B15) 
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